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CHAPTER  I 


INTRODUCTION 


i*  Historical  background 

Although  the  history  of  the  kinetic  theory  of  gases  lias  been  traced 

1* 

back  to  ancient  Greece,  its  scientific  development  started  with  other 

branches  of  modern  science.  The  first  main  success  came  in  1859  when 

2 

Mucwell  discovered  the  law  of  the  distribution  of  the  molecular  veloci¬ 
ties  for  a  gas  in  equilibrium  and  rediscovered  the  equipartitlon  of  the 
mean  molecular  energy  and  as  a  consequence  the  ideal  gas  laws. 

Definite  advances  were  made  in  1879  by  Boltzmann^  who  discussed  es¬ 
pecially  the  approach  to  equilibrium.  He  established  the  famous  integro- 
differential  equation  (the  Boltzmann  equation'  which  the  one  particle  dis¬ 
tribution  function  must  satisfy  whatever  the  state  of  the  gas  is.  There¬ 
after  one  of  the  main  problems  has  been  how  to  solve  this  equation  and 
how  to  derive  macroscopic  equations  such  as  the  hydrodynamlcal  equations. 

In  1917  Enskog  published  his  Uppsala  dissertation  in  which  he  gave 
a  general  method  for  the  determination  of  the  distribution  function  from 
the  Boltzmann  equation.  His  method  was  a  modification  of  a  method  first 
proposed  by  Hilbert.-*  He  derived  the  general  formulae  for  the  viscosity, 
heat  conduction  and  diffusion  of  gases.  At  almost  the  same  time  Chap- 
man^'^  obtained  independently  the  identical  results  with  a  slightly  dif¬ 
ferent  method. 

•Raised  numbers  refer  to  the  bibliography  on  pp.  96  and  97. 
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A  very  important  extension  of  the  theory  was  made  by  Enskog'*  in 
1922.  Using  the  elastic  sphere  model  for  the  molecular  interactions,  he 
p  w^ed  how  the  theory  could  be  modified  for  dense  gases.  If  one  regards 
the  Boltzmann  equation  as  the  description  of  the  state  of  an  ideal  gas 
which  is  not  in  equilibrium,  then  the  Enskog  theory  may  be  considered  as 
an  attempt  to  describe  the  non -equilibrium  properties  of  a  non-ideal  gas, 
according  to  the  ideas  of  van  der  Waals  (1873)* 

During  the  later  part  of  these  developments  of  the  non-uniform  gas 
theory,  the  statistical  mechanics  of  the  equilibrium  state  was  firmly  es- 
tablished  on  the  ensemble  theory  by  Gibbs.7  In  1937#  Mayer  succeeded 
in  deriving  the  equation  of  state  for  dense  gases  with  arbitrary  central 
molecular  forces  from  the  point  of  view  of  statistical  mechanics.  It  is 
therefore  understandable  that  in  the  next  decade  the  general  trend  of  the 
kinetic  theory  was  the  elucidation  of  the  connection  of  the  Boltzmann 
equation  and  of  the  Enskog  theory  with  the  Llouvllle  equation,  which  is 

the  basis  of  the  ensemble  theory.  The  first  investigations  in  this  di- 

^  11 
rection  were  dene  mainly  by  Born  and  Green,  and  by  Kirkwood  and  his 

Ip 

collaborators.  They  introduced  the  higher  order  distribution  functions, 
and  found  the  hierarchy  of  equations,  which  will  be  discussed  in  Chapter  II. 

2.  Assumptions  in  the  Boltzmann  Equation 

Before  discussing  the  kinetic  theory  further,  it  is  important  to 
know  the  basic  assumptions  on  which  the  Boltzmann  equation  depends.  It 
is  well  understood  now,  that  there  are  two  such  assumptions.  The  first 
one  is  the  assumption  of  binary  encounters.  Since  the  molecules  of  a  gas 
interact  with  each  other  through  a  short  range  force  (range  r0  is  of  the 
order  of  10*®  cm),  for  a  dilute  gas  the  average  volume  per  perticle  will 
be  large  compared  to  the  volimrn  of  the  action  sphere  ksT^ 3.  Therefore 


3 


the  probability  of  finding  two  molecules  in  the  same  action  sphere  is  very 
small,  and  the  probability  of  finding  three  or  more  molecules  interacting 
simultaneously  will  be  quite  negligible.  The  neglect  of  such  triple  in¬ 
teractions  amounts  to  the  assumption  of  binary  encounters .  Clearly  for 

dense  gases  this  assumption  will  have  to  be  modified. 

«• 

The  second  assumption  is  the  so-called  ’’Stosszahlansatz*'-3^  about  the 
number  of  pairs  of  molecules  vhich  are  in  the  position  to  collide  during 
a  given  short  time  interval.  It  is  well  known  that  this  statistical  as¬ 
sumption  makes  the  equation  irreversible  in  time.  Since  the  gas  consid¬ 
ered  as  a  mechanical  system  is  reversible  in  time,  it  is  clear  that  the 
time  used  in  the  Boltzmann  equation  is  not  the  exact  mechanical  time,  but 
is  measured  on  a  coarser  scale,  in  which  in  each  time  element  a  great  num¬ 
ber  of  collisions  occur.  While  the  Stosszahlansatz  seems  very  plausible 
and  is  verified  by  many  true  consequences  of  the  Boltzmann  equation,  it 
is  desirable  to  replace  the  Ansatz  by  more  general  statistical  assumptions, 
so  that  the  extension  to  triple  and  higher  orfcer  collisions  would  become 
possible . 

In  the  attest  to  derive  the  Boltzmann  equation  from  the  Liouville 
12a  »b 

equation,  Kirkwood  had  to  average  the  latter  equation  over  a  very 
short  time  of  the  order  of  the  collision  time.  In  his  theory,  only  the 
binary  collisions  are  considered  and  the  Stosszahlansatz  is  replaced  by 
a  new  assumption.  Kirkwood  as  suae  s  that  if  the  two  particles  involved  In 
a  binary  encounter  are  far  apart  from  each  other,  the  binary  distribution 
function  is  the  product  of  the  one  particle  distribution  functions  at  the 
corresponding  positions.  This  assumption  seems  plausible  when  the  gas  is 
not  very  far  from  the  state  of  the  local  equilibrium  and  it  turns  out  to 
be  equivalent  to  the  ’’Stosszahlansatz"  when  the  spatial  distribution  of 
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the  molecules  of  the  gas  varies  very  slowly.  For  denser  gases,  where 
triple  collisions  become  important,  we  shall  see  that  it  will  be  neces¬ 
sary  to  modify  the  Kirkwood  assumption. 

14 

3 .  The  general  idea  of  Bogolubov 

At  the  same  time  with  Born  and  Green  and  Kirkwood,  Bogolubov  pro¬ 
posed  a  more  satisfactory  theory.  This  theory  can  be  interpreted  in 
various  ways  (see  for  instance,  reference  15);  we  will  try  in  the  follow¬ 
ing  to  show  its  relation  to  the  Kirkwood  idea  of  time  averaging  or  "coarse 
graining  in  time." 

There  are  three  features  which  are  characteristic  for  the  Bogolubov 
theory.  The  first  one  is  the  use  of  successive  time  scales  of  increasing 
roughness.  Then,  it  will  always  be  assumed  that  the  theory  describes  the 
state  of  a  gas,  which  does  not  deviate  very  far  from  a  local  equilibrium 
state.  And  finally  all  properties  of  the  gas  will  be  expanded  in  powers 
of  the  average  concentration,  analogous  to  the  virial  development  used  in 
the  theory  of  the  equilibrium  state.  In  this  way,  triple  and  higher  order 
collisions  are  successively  taken  into  account. 

Imagine  at  time  t  *  0  a  severely  disturbed  state  of  the  gas  very  far 
removed  from  the  equilibrium  state.  The  temporal  development  of  the  state 
of  the  gas  could  then  be  described  only  by  the  Liouville  equation.  We  as¬ 
sume  that  after  a  very  short  time  of  the  order  of  the  collision  time 
•12 

t#**10  sec,  the  state  of  the  gas  relaxes  to  a  quasi -equilibrium  or 
"normal"  state,  in  which  the  description  of  the  state  of  the  gas  can  be 
simplified.  The  effects  of  the  intermolecular  forces  on  the  one  particle 
distribution  function  are  c’toothed  out,  and  hence  it  will  vary  little  in 
a  time  of  order  rQ.  But  the  binary  and  higher  order  distribution  func¬ 
tions  still  suffer  the  direct  effect  of  the  intermolecular  forces  and 
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will  change  rapidly.  However  if  we  average  these  functions  successively 
over  times  of  order  tc,  the  remaining  time  dependence  of  these  functions 
would  be  due  to  the  change  of  the  one  particle  distribution  function. 
Therefore  on  the  first  "coarse  grained"  time  scale,  which  we  will  call 
the  kinetic  time,  the  higher  order  distribution  functions  would  depend 
on  time  only  through  the  one  particle  distribution  function. 

In  this  "kinetic"  stage,  all  the  higher  distribution  functions  will 
be  expanded  in  powers  of  the  average  concentration.  Actually  the  expan¬ 
sion  parameter  i6  the  average  density  of  particle  .times  r*,  which  will  be 
small  if  the  gas  is  far  from  the  condensation  point.  One  then  assumes 
that  the  lowest  order  term  is  the  product  of  the  one  particle  distribu¬ 
tion  functions  in  certain  phases  of  the  states  of  the  particles.  This  is 
the  modification  of  the  Kirkwood  assumption  which  can  be  used  for  denser 
gases . 

After  a  time  of  the  order  of  the  time  between  collisions  tQ  At  10"8 
sec,  even  the  description  using  the  kinetic  time  scale  becomes  unneces¬ 
sarily  detailed  for  mo3t  purposes.  There  are  exceptions;  for  instance  for 
strong  shockwaves,  the  kinetic  time  description  is  required  to  investi¬ 
gate  the  almost  discontinuous  changes  of  the  properties  of  the  gas.  How¬ 
ever  in  most  cases,  the  gas  is  so  close  to  local  equilibrium  everywhere 
that  the  change  of  the  macroscopic  quantities  (density,  temperature  and 
the  macroscopic  mass  velocities)  is  slow  coaq>ared  to  the  detailed  change 
of  the  first  distribution  function.  Therefore  if  one  averages  successive¬ 
ly  by  over  times  of  order  t0,  then  on  this  second  coarse  grained  time  scale 
(which  we  will  call  hydrodynamic  time ) ,  one  can  say  that  the  one  particle 
distribution  function  depends  on  the  time  only  through  the  macroscopic 
quantities. 
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The  development  of  the  theory  in  this  "hydrodynamic"  stage,  turns 
out  to  be  very  similar  to  the  theory  of  Chapman-Enskog .  All  quantities 
are  again  expanded  in  powers  of  a  parameter  p,  which  is  a  measure  of  the 
uniformity  of  the  macroscopic  quantities.  Actually  p  is  of  the  order  of 
the  relative  change  of  these  quantities  over  a  mean  free  path. 

♦ 

k.  The  purpose  of  the  present  dissertation 

The  purpose  of  this  dissertation  is  to  elucidate  the  Bogolubov  the¬ 
ory  as  much  as  possible  and  to  carry  out  the  theory  till  the  hydrodynami- 
cal  stage  (which  Bogolubov  only  indicated),  in  order  to  obtain  the  formal 
expressions  for  the  viscosity  coefficients  and  heat  conductivities  of  a 
dense  gas.  There  are  two  main  problems.  The  first  one  is  the  solution 
of  the  equations  of  motion  for  three  or  more  interacting  particles.  The 
second  problem  is  to  obtain  the  distribution  functions  assuming  that  the 
first  problem  is  solved.  We  will  be  concerned  only  with  the  second  prob¬ 
lem.  It  turns  out  that  the  two  problems  can  formally  be  separated,  and 
that  it  is  possible  to  find  the  deviations-  of  the  distribution  functions 
from  their  local  equilibria  forms  in  terms  of  the  formal  solutions  of 
the  particle  mechanics. 

jj.  The  results  obtained 

1)  The  relation  between  the  Bogolubov  theory  and  the  Kirkwood  theory 
is  clarified. 

2)  The  effect  of  triple  collisions  of  particle  is  estimated  in  the 
kinetic  theory. 

3)  The  theory  of  the  hydrodynamical  stage  has  been  developed  up  to 
the  second  order  in  the  uniformity  parameter  p. 
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a)  In  the  first  order  of  the  uniformity  parameter,  one  obtains  the 
ideal  or  Euler  hydrodynamical  equations,  in  which  the  pressure  as  func¬ 
tion  of  the  density  has  the  same  form  as  in  equilibrium  according  to  the 
Mayer  theory. ^  /Iso  the  energy  equation  is  derived,  which  in  this  ap¬ 
proximation  corresponds  to  adiabatic  changes,  and  in  which  the  energy 
density  is  again  the  same  as  in  the  equilibrium  theory. 

b)  In  the  second  order  of  the  uniformity  parameter,  one  obtains  the 
Stokes -Navier  equations.  There  occur  two  viscosity  coefficients,  for 
which  one  obtains  expansions  in  the  number  density  of  the  forms: 


(0)  (1)  2  (2) 

qx  -  m  +  n  rji  ♦  tjx  + 

(1)  2  (2) 

Tj2«  n  +  r.  rj2  + 


(1.1a) 

(1.1b) 


where  the  Tji  is  the  Chapman -Enskog  value  for  the  rhea/  viscosity  co¬ 
efficient,  and  1)2  is  the  bulk  viscosity.  In  the  energy  equation  a  term 
appears  corresponding  to  heat  conduction,  and  for  the  heat  conductivity 
one  obtains  the  analogous  expansion: 


(0) 


♦  n  t 


(1) 


+  n 


(2) 


(1.2) 


where  is  again  the  C hapman -Ens kog  value  for  the  heat  conductivity 

coefficient. 

4)  *When  the  molecules  of  gas  are  hard  spheres,  the  parts  of  the 
formal  coefficients  and  which  depend  on  the  binary  collisions 

can  be  evaluated,  and  these  parts  agree  up  to  the  first  order  in  n  with 
the  Enskog  theory  of  dense  gases.  Even  for  this  simple  model  the  complete 
first  order  terms  h-'ve  not  been  found  because  of  the  difficulty  of  the  me¬ 
chanics  of  triple  collisions. 


CHAPTER  II 


GENERAL  DESCRIPTION  OF  THE  SYSTEM 

1.  The  Llouvllle  equation 

In  the  following,  we  shall  consider  the  behavior  of  a  system  of  N 
identical  molecules  in  a  vessel  (volume  V),  which  obey  the  law  of  clas¬ 
sical  mechanics.  To  simplify  the  problem,  we  restrict  ourselves  to  point 

molecules  repelling  each  other  by  a  known  monotonic  central  force  poten- 

* 

tial  0  between  each  pair  (i,J),  which  is  a  function  of  the  distance  be¬ 
tween  the  pair  only  and  which  has  a  finite  very  small  range  rQ,  so  that 
0(0)  ♦  »  and  0(ro)  *  0.  Furthermore,  the  system  is  supposed  to  be  not 
under  any  outside  force  except  the  force  due  to  the  wall  potential  of 
the  container. 

Let  the  coordinates  and  momenta  of  the  i-th  particle  be  *  (<lj>Pi)- 
The  state  of  ".he  system  at  time  t  is  completely  determined  by  the  set  x^, 
x2  ...  xN.  Usually  it  is  convenient  to  introduce  the  6N-dimensional  phase 
space  for  the  system  as  a  whole,  the  T -space. ^  The  state  of  the  system 
is  then  represented  by  a  point  In  this  space,  and  the  temporal  develop¬ 
ment  of  the  system  is  completely  represented  by  the  trajectory  of  this 
point . 

For  a  system  with  a  large  number  of  particles  it  is  physically  not 
meaningful  to  assign  the  initial  state  completely,  since  only  some  average 
values  corresponding  to  the  results  of  macroscopic  measurements  are  known. 
Therefore  one  must  consider,  in  the  language  of  Gibbs ?  an  ensemble  of 
identical  systems  differing  in  their  initial  states  and  follow  the  streaa- 
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ing  of  the  "ensemble  fluid"  in  time.  In  other  words,  one  must  consider 
a  probability  distribution  . .  .xN;  t)  in  the  T-space  and  follow  its 
development  in  time. 

According  to  the  definition  of  the  probability  distribution  function. 


/-/ 


Dn(x1,...xn;  t)  dxj,. .  ,dx 


N 


(2.1) 


using  the  notation  dxA  *  dq^dpi .  The  change  in  the  probability  distribu¬ 
tion  function  with  time  is  determined  by  Liouville’s  equation: 


-  {%.  Df} 


(2.2) 


where  HH  is  the  Hamilton  function.  For  the  system  considered,  it  is  gi¬ 
ven  by: 


Hn  * 


(2.3) 


where  m  is  the  mass  of  each  molecule,  •  0(  I )  the  interaction 

potential,  and  ^(q*)  is  the  potential  produced  by  the  walls  of  the  vessel, 
so  uhat 

.h  inside  the  vessel, 
walls  of  the  vessel. 

The  Poisson  bracket  DjjJ  can  be  written  as: 


M?i) 


r  o  if  1’ 

[+*  at  the  > 


{V"«] 


f  _  an,  ao, 

1“1  dp^  dl>i  5q 


t) 


(2.4) 


where  for  any  vector  |  we  have  written  b/b |  «  grad|  .  Introducing  the 
"Hamiltonian"  operator 

s  ♦  s 


*•<* .  I  h.  z 


1*1  1  i<J 


(2.3) 
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with 


%j  •  .  J_  ♦  $4 . 4- 

&ti  d£i  dpj 

one  can  put  the  Liouville  equation  (2.2)  in  the  forms 


dPy 

at 


(2.6) 


(2.7) 


if  the  particles  are  inside  the  vessel. 

The  Liouville  equation  is  an  immediate  consequence  of  the  equations 
of  motion  and  allows  us  in  principle  to  find  Djg(xx . .  .x^;  t)  if  the  initial 
distribution  Djj(xx . . .xjj;0)  is  given.  Since  all  particles  are  indistin¬ 
guishable,  one  must  choose  Djj(xx . .  .xN;0)  as  a  symmetric  function  of  xx . . . 
Xfl.  For  this  initial  distribution,  Djj(xx . .  .Xjjj  t)  remains  symmetric  be¬ 
cause  Hjj  or  ^(Xx..  .Xjj)  are  symmetric  with  respect  to  all  particles. 


2 .  The  B-B-G-K-Y  equations 

For  very  large  N,  it  is  practically  impossible  to  obtain  an  explicit 
expression  for  Djy(xx . .  .xj|f  t),  because  this  involves  the  precise  integra¬ 
tion  of  motions .  This  is  also  not  required  since  one  is  only  Interested 
in  the  change  of  some  macroscopic  quantities  with  time. 

The  macroscopic  quantities  which  have  the  most  direct  physical  mean¬ 
ing  are  for  a  small  volume  around  some  point  in  spacer  the  number  of  particles, 
the  average  velocity  of  the  group  of  particles,  the  average  total  kinetic  en¬ 
ergy  cf  the  group  of  particles,  and  the  average  total  energy  of  the  group  of 
particles.  Fortunately  these  quantities  depend  not  an  the  complete  distribu¬ 
tion  function  Dg(xx . .  .xj|{  t)  in  T-space,  but  on  the  probabilities  of  find¬ 
ing  a  single  particle  in  a  certain  range  dq  dp  around  a  phase  point  (q,p), 
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or  of  finding  two  arbitrary  particles  in  certain  ranges  dqx  dpx  and  dq^  dp2 
around  two  phase  points  (qi,px)  and  (q^iPa)  irrespective  of  the  phases  of 
all  other  particles.  These  probability  functions  are  obtained  by  integra¬ 
ting  Dj|(xi . .  .xj^;  t)  over  all  except  one  or  two.  Since  D^(xx...xn;  t) 
is  symmetric  in  xx,...,Xtf,  these  functions  will  be  the  same  whatever  single 
or  pair  of  particles  is  chosen. 

By  integrating  the  Llouville  equation,  one  obtains  a  hierarchy  of 

14 

equations  derived  independently  and  simultaneously  by  Bogolubov,  Born 
11a  12  17 

and  Green,  Kirkwood,  and  Yvon  end  therefore  called  the  B-B-G-K-Y 
equations.  Following  Bogolubov,  introduce  the  partial  distribution  func¬ 
tions  by: 


^  F8(x1,...x8jt)  «  J'' 4 'J  %(xi  •  •  •Xjj»t)dx8+1. .  .dxjj  (2.8) 

In  the  limit  U  ♦  »,  V  ♦  »  end  v  a  v/N  finite,  if  one  considers  only  those 
molecules  deep  inside  the  vessel,  one  may  forget  the  effect  of  the  walls. 
Therefore,  integrating  Llouville 's  equation  in  the  form  given  by  (2.7) 
over  *a+l >  *  * ' >XN  and  multiplying  by  V8,  one  obtains  immediately: 


&Fg 

*5 r 


♦^FS  *  V*]''  *  ^*6+1 '  •  ,dxH  ^  —  -S“ 

L  i«S+l  ®  «$i 


Z 

•♦1 <XN 


hi  * 


selOeCI  J 

mm 


Since  D^(xx,...xg;t)  must  be  assumed  to  vanish  for  large  |pi|  and  |f^| 
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I  •••/  d*B.l"-dxN  ®kl  •**!*)  -  0 


1  >  s+1 


k  >1  >  S4l 


For  the  term: 


V°  l<s  I  "'I dX8+1" -,iXli  9lJ  D»(Xl •  •  •x«'t> 

B+iijSn 

all  the  contributions  from  different  values  of  J  are  the  same  because  of 
the  symnetry  of  Djj(xi,  . .  .xjgjt ).  Therefore  using  the  definition  (2.8) 
again: 

Vs(N-S)  L  J dxs+1  fti84  J dxs+2-..<ixR  DN(xi...xN;t)  » 

N-s  y  i 

“  T"  i <SJ  <^S4l®is4lF84l(xl-***B4lit) 


which  in  the  limit  stated  above  becomes : 


7  i<B  J^B+l  ®iS4l  F84l  (*l***XS+l»i') 


for  fixed  s.  Hence: 


^  •  7  J^B.l  ®1.+1  FB.l  8  •  l>2> 


(2.9) 


Later,  we  shall  be  especially  Interested  in  the  cases  s«l  and  s*2  which 
become,  using  the  explicit  forms  of  Xi  sod  Xe. 

frilgiiM  4  jj^FifcxJt)  .  i  f  to*  *ia  F*(xlfXatt)  (2.10) 

•  £  J Fs(*i »*•**»!*)  (2.11) 


CHAPTER  III 


THE  BOGOLUBOV  THEORY  OF  THE  KINETIC  STAGE 

1.  The  b&alc  equations 

Ab  explained  In  the  introduction,  we  will  assume  that  on  the  first 
coarse  grained  time  scale  (the  kinetic  time),  the  higher  order  distribu¬ 
tion  functions  depend  on  time  only  through  the  first  distribution  func¬ 
tion  Fx,  so  that: 

Fs(x1...xs,t)  - *  F8(xx...x8|Fx)  (A) 

where  the  vertical  bar  denotes  that  Fs  depends  functionally  on  Fx.  The 
whole  time  dependence  sits  in  Fx,  and  the  form  (A)  is  assumed  to  be  valid 
for  any  Initial  distribution  D)g(xx. . .x^jO)  after  an  initial  period  of  or¬ 
der  t0.  The  first  distribution  function  is  expected  to  vary  smoothly  on 
the  kinetic  time  scale  and  to  fulfill  the  basic  kinetic  equation  of  the 
form: 

•  A(X»|FX)  (B) 

St 

The  unknown  functionals  F8(xx. . .x8 |FX)  and  A(xx|Fx)  must  follow  from  the 
hierarchy  of  equations  (2.9) •  To  determine  them  in  successive  approxima¬ 
tion,  we  develop  both  functionals  in  powers  of  l/v  (vlrial  expansion): 

+M,<*ip+£A*fciir>+-  -  (j.i) 

E <*<  •  -  Iv<V  %ff?) + .  (j.*) 


V 
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Ik 


Substituting  the  series  expansion  for  F2  in  the  first  equation  (2.10) 
of  the  hierarchy  one  obtains: 


v'j&ft.F/ 


and  comparing  this  with  the  kinetic  equation  (B),  ve  get  immediately: 


(3-3a) 


and  in  general 

At(*‘ 1 P  -  pt  8,k  x,  i  p )  (5.3c) 


Since  Fs(xx. . .x8 |FX)  depends  on  tine  only  through  Fx,  one  can  egress 
its  change  in  tine  by  the  kinetic  equation  (B).  Let  the  first  order  var¬ 
iation  of  F8(xx.  .  .x8  |FX)  for  the  variation  o*  Fx  to  Fx  +  c6Fx  be 
(xx...x8  |FX,6FX),  then  obviously,  ^8(xx.  *  .*s  lFi»6Fx)  18  linear  in  the 
(6Fx)'s  which  will  have  different  arguments.  The  argument  X  of  6Fx(X{t) 
is  determined  by  the  functional  form  of  Fs(xx. . .x8|Fx).  Replacing 
6Fx(X;t )  by  A(x|Fx(jt))  one  gets  the  change  of  Fs(xx. .  .x8  |FX)  in 
time.  Writing: 


(3.W) 


On  substituting  the  l/v  expansions  of  Ft(xx...xs  |FX)  and  A(|FX) 
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'H 


(3-i.b) 


where  the  operators  acting  on  any  functional  t(xx..>  |fx ) 
fined  by: 


of  Fx  are  de- 


(5-5) 


Comparing  (3*^)  with  the  B-B-G-K-Y  equations  for  s  >  2,  using  the  series 
expansion  for  Ffl+^  and  equating  equal  powers  of  l/v,  one  obtains: 


*1 


(3.6a) 


(3.6b) 


(3.6c) 


These  equations  are  functional  equations  for  Fx,  so  they  must  hold  for 
any  function  Fx. 

jj.  Determination  of  the  functionals  FB*^(xx...xB |FX) 

Tie  program  for  solving  the  basic  equations  (3*3a,to,c)  and  (3.6a,b,c ) 

is  as  follows.  The  first  equation  (3.6a)  should  determine  ?i°^ i  this  In 

turn  determines  Ax  according  to  (3-36),  so  that  also  DXFB°^  Is  known. 

Hence  in  (3>6b)  the  right  hand  side  is  known,  and  (3.6b)  should  determine 
(1) 

F,  ’ ,  from  which  A*  follows,  and  so  cm.  In  th?s  section  we  will  there¬ 
fore  assume  that  the  right  hand  sides  of  the  equations  (3.6e,b,c)  are 
known,  and  ask  how  the  unknown  functionals  FB*  ^  can  be  found. 
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It  should  he  emphasized  that  equations  (3.6a,b,c)  contain  functional 

derivatives.  Therefore,  just  as  with  first  order  differential  equations, 

sane  sort  of  initial  or  boundary  conditions  arc  needed  In  order  to  deter* 

mine  the  solutions  completely.  To  formulate  these  conditions,  we  Intro* 

(k) 

duce  the  streaming  operator  '(xx. . .xk)  by: 

(*,  **)  -  e*f>{TXK<r'-r,)J  (5’7) 

It  Is  the  time  displacement  operator  (over  time  t)  in  the  streaming  of 
the  k  particles  In  their  phase  space  under  the  influence  of  their  mutual 
interactions.  One  easily  sees  that  If  (q±  p*)  Is  the  phase  of  particle 
1,  then  in  the  motion  of  the  k  particles  the  phase  of  particle  1  a  time 
r  later  is  sWfil-  Also  when  X(xj,.. .  .xjJ  Is  an  arbitrary 

function  of  the  phases  of  the  k  particles,  then: 


(3.8) 


Clearly  the  sjk)  forms  an  additive  Abelian,  one  parameter  group  of  op¬ 
erators  with  parameter  t;  4J’  •  -  4iir* »  sij’  is  the  inverse 

of  4k>  and  »  1. 

We  will  now  assume  following  Bogolubov  that  for  any  of  the  function¬ 
als  F«: 

s-a.j/—  (c> 

The  basic  nature  of  this  ass  taction  should  be  e^hasized.  It  replaces 
and  generalizes  the  Boltzmann  8 toe szahlansstz  and  the  Kirkwood  assumption. 
It  expresses  somehow  the  requirement  that  outside  some  action  volume  the 
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correlation  between  the  particles  vanishes.  This  Is  especially  clear 
when  the  system  is  spatially  uniform.  Since  In  general: 

-  Fcftlr,  f;t) 

clearly  In  the  spatially  uniform  case  S^Fxfxjt)  »  Fx(xjt).  Because 
there  Is  a  strong  repulsion  between  the  particles,  the  relative  distances 
of  the  s  particles  In  the  phase  S^Xx,  . ..  si?^x8  will  be  large  If  In 
the  phase  Xx*  . . .  xB  they  were  in  each  others  action  spheres.  One  can 
also  say  that  the  two  basic  assumptions  (A)  and  (C)  for  the  s  particle 
distribution  function  are  the  two  properties  of  the  equilibrium  distri¬ 
bution,  which  one  assumes  that  F8  already  has  in  the  kinetic  stage.  In 
Section  6  we  will  actually  see  that  (C)  is  fulfilled  in  equilibrium. 

Using  the  virial  development  for  F8,  (C)  implies  that: 

(CD 

and 

^  IS? F)  -0  l  *0  <C*) 

TO  find,  with  the  condition  (Cl)  f£°),  one  replaces  in  (3.6a)  the 
functional  derivative  by  an  ordinary  derivative  in  the  following  way. 

Put,  in  D0F^°)(xx...*8|Fi),  the  function  S^Fx  for  TXf  then  by  the 
definition  (3-9)  of  the  operator  D*, 

Suppose  the  variable  of  one  of  the  |sj1^F1)  in  the  bracket  is  X  a  (?,?), 


then: 
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A.<XI$F>  - 
.  - 


n>r 

■>di 


-Rot*) 


by  (3-3a) 


by  (3.7) 


-  -2£t 

*>r 

since  Fx(Xjt)  is  Independent  of  t.  From  (3.1*a): 

r  Sgval^Q  .  'X+t  l  -  ^EV-r.ist’F) 

L  S(S^)  rr  J  ->r 

Hence,  considering  the  linearity  of  the  functional  derivative: 

(■RTTV-at.lf;)], 


*  - 

'5rF-F 


On  replacing  Fx  appearing  in  (3.6a)  by  sj^'Fi, 

>r 


^  -  Ti^  vT^r,-  *sl  £'ff )  -  O  (5-9) 


then,  the  "solution"  of  the  equation  is: 


(3.10) 


because  of  (3.7). 


Operating  with  the 


from  the  left  gives: 


This  is  the  condition  which  the  functional  fox*  of  F^0^(xx...x>|F|>) 
should  satisfy.  In  the  above  equation,  the  left  hand  side  is  indepen¬ 
dent  of  t,  so  this  aust  hold  for  an  arbitrary  r.  Therefore,  taking 

T  ♦ 
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£V-*i1E  )  *  $.r 


-  JL  sX 

'T-*e» 


(3.11) 


according  to  (C£). 

To  find  next  fP^Xj..  .  .xs  |FX),  put: 


•  (3.12) 


then  (3*6b)  becomes: 


id  -  P*’  -z‘iF,)  <5-«> 


After  the  same  calculation  as  before: 


faFfo-siF)]  -  - 

1  r*  r’ Js?f>F?  *»r 

and  from  the  functional  equation  (3*13)  For  Fx,  (replacing  Fx  by  s£^Fx): 

2ghpl£El  _  #tEV*w?fj (j.u) 


Putting 


ffi  -x.lSpfr  ;?■!£) 


(3.15a) 


(3.15b) 
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then  clearly: 


(3.16) 


and  one  finds: 


Q<ff3r>-3;;r Ip 

>r 


-&t 


This  may  be  integrated  easily,  yielding: 


from  which  follows: 


j^a-,"-*slK)  *  l2£K) 

Again  the  left  hand  side  is  independent  of  t,  therefore  taking  the  limit 


t 


Ff(vnlF)  -  £;*?"£(*.••• ■*«?& >+prS!?#«~',lS?F> 


fir  Si?  fo'.-z.  |S?K) 


(3-18) 


according  to  (C$). 

The  higher  fJ^  are  found  in  the  same  way.  One  gets: 

CV*.lE>  » pr£"  #V-*,|£jv> 


where: 


(3.20) 
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1‘  ^termination  of  the  functionals  Ai(x|Fx) 

According  to  (3.5b)  Ax(xxjFx)  Follows  from  F^(xxxa |FX),  which 
we  write  in  the  form: 


(5.21) 


where  the  subscripts  +®  Imply  the  limit  indicated  in  equation  (3.11), 
and  where  for  abbreviation: 


$ 


\<h~',  I) 


(3.22) 


To  separate  the  mechanics  from  the  distribution  functions,  it  is  often 
convenient  to  write  instead  of  (3.21): 

K  (3-23) 

where  corresponding  t'  x*  •  (q^Pi),  Ci  •  (ti,ni)  and  5  is  the  Dirac 
6-function.  From  (3. 3b)  one  then  can  write: 

mjjTt  9„ f?(i, ■>*>£(*, ■■*) 

(3—  *) 

with 

Q/r.  1  <■<,'  -  fa  f  $<*-<-  )  (5.25) 

i*t 

Turning  now  to  Aa(xx|Fx),  we  have  first  to  find  Fj^(xxxa|Fx). 


According  to  the  definition  of  the  operator  Dx: 
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=  Sl+t.  (  JjVrtoAfaip+AftJfrFiwt)] 

-  S5rO>p*.  [  iS%’‘S>>+a,2%,>s&j  jfa  :t) 

where  we  have  omitted  x's  in  the  arguments  of  the  operators  s£^. 
according  to  (3.11); 

F'V.^lf’)  *  S-L( r< ** r* )  ft,  Z»}F( *i ; t) 

*#  fc*» 

-  SUL"  * TTK'ir.i) 


Since 


one  finds: 


frvJf-)  - 


•ini)  - 


~  » +^j  Si''  »<&'»)]  7T  fat) 


(3.26) 


Using  equation  (3*18)  one  then  can  write  F^fx^XalF*)  in  the  form  an¬ 
alogous  to  (3.23): 

£  ^ J rr4  •  >  (3.27) 

where  [  ]  is  the  operator  occurring  between  square  bracket  in  (3.26). 

Note  that  this  is  syoaetrlc  in  Xj,  and  xg.  Finally  from  (3.3c): 

A(,*iFT>  'jjfj^AJs,  0.28) 


with 
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rV’-KW  J  gjfaf  kri-to 


(3.29) 


k .  The  "size"  of  Aa(xi|Fi) 

The  ratio  of  A2  to  Ax  should  be  of  order  r§,  if  rQ  is  the  range 
of  the  interaction  potential,  so  that  the  development  parameter  in  the 
kinetic  equation  is  really  rj*/v.  This  will  be  so  for  dimensional  reasons 
provided  that  the  integrals  in  Aa  are  convergent  and  extend  only  over 
the  action  volume  of  the  three  particles  xx,  xg,  x3. 

To  see  this,  consider  ft3(xx |5x£aC3) .  The  integrand  of  the  x2  in¬ 
tegration  will  vanish  if  |q2-qx |  >  rQ  because  of  the  A12  operator. 
Consequently  one  needs  to  consider  only  those  x2  for  which  jqa-qxj  < 
r0.  For  this  phase  (xx,x2)  the  S^(xxx2)  operator  separates  the 
.  particles  1  and  2  by  a  distance  r0  in  the  time  r  &  rQ,  and  from  that 
time  on  the  distance  of  the  particles  Increases.  Putting: 

■($.£)  t  -  i.  a. 

let  r*  be  the  time  for  which  |^a-qx  |  %  2r0.  Now  the  operator  occurring 
in  the  t -integral  can  be  written  as: 

jdT,  S.'li*  *.  *  v  - 

-  li* (5.30) 

For  t  >  0x*(Xx2*)  *  0  and  for  fixed  X3  at  least  one  of  the 

•ts(^i»*s)  “4  %a(^e>*  )  suet  be  tero,  since  xx  and  x*  sure  sepa¬ 
rated  by  a  distance  bigpr  than  2r0.  S\gppose  ®xj(xi,x3)  -  0,  but 


2k 


©23 $1  **3)  f  0.  Prom  the  definition  of  S  ^(*1*2X3),  one  then  has: 


S&UV  —  Stfrv&i.-. 


and  since  in  this  case  6<?><8A)  s‘la)(SA>  *  1  one  sees  that  the  op¬ 
erator  (3*30)  goes  to  zero.  This  is  also  the  case  if  ©saC^Xg)  *  0 
but  ©is^ifXs)  /  0,  so  that  one  can  conclude  that  the  t- integral  will  be 
convergent,  and  that  the  integrand  will  only  be  different  from  zero  for 
a  time  of  the  order  of  a  collision  time. 


2*  Spatially  uniform  systems 

When  Fitxjjt)  does  not  depend  orrthe  <J1#  i.e.,  when  the  system  is 
spatially  uniform,  it  is  possible  to  simplify  the  expressions  for  Ax(xx|Fx) 
and  AgfxxlFx). 

We  will  show  that  in  this  case  Ax  can  be  reduced  to  the  Boltzmann 
collision  integral. 

Proof:  Since  F*  is  independent  of  the  spatial  coordinates  Fx(CxJt)  * 
a**4  *  Fx(Ti2jt).  In  equation  (3.210  one  can  then  per¬ 
form  the  l-integratlons,  and  since  the  operator  does  not  change  the 

momenta,  one  gets: 


with: 

$(/H, )$$•-*),)  (}•»> 

Let  S^^(*x**)?i^i2^  i  •  1»2  then  it  is  easily  proved  that  the 
are  functions  of  $x»  $e  and  of  the  relative  coordinate  r«x  •  q*-qx • 
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Since  ©j,2  also  depends  only  on  the  relative  coordinates,  one  sees  that 
by  taking  in  the  X2-integral  the  origin  for  ^  in  that  n2  and  there¬ 
fore  Ai  will  depend  only  on  jpj.,  as  it  should.  We  will  write  instead 
of  (5.32): 

iU>  -\K\ dft  6L5*f^>  -1,v H'  1) 

Kj^«  ft.  )  (3.33) 

*(2) 

The  are  the  constant  initial  impulses  in  the  binary  collision 

governed  by  the  Hamiltonian  Hg  which  leads  to  the  phases  xi,x2  at  time 
zero .  Therefore : 

H.  -  t')+L  -  iir  f f  fH + ( f*T) 

Hence  according  to  the  definition  of  the  Poisson  bracket: 

{  h*(  -  0 

or 


In  (3*33)#  take  for  the  r^  integration  cylindrical  coordinates  with  the 

axis  in  the  direction  of  the  relative  velocity  g  •  i  (p2-Pi).  Call  the 

■ 

coordinate  along  this  axis  I,  and  the  polar  coordinates  perpendicular  to 
the  axis  (#,0).  Then: 
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where  g  *  |g|.  The  i-integration  can  be  done  Immediately: 


f 

scff-fyscfy-l) 


At  i  »  -«o  the  two  particles  are  outside  their  interaction  region,  and 
since  the  s(2)(xxx2)  operation  will  never  bring  them  together 


JL  p 
$*-00  1  ‘ 


(*) 


ft 


i  «/,  I 

On  the  other  hand  for  t  *  +»,  although  the  particles  are  then  also  out¬ 
side  their  interaction  range,  the  s(£)(xjx2)  operation  will  produce  a 
collision.  Therefore: 


fe.  ?'*'  -  f:  i-..* 

where  the  are  the  iapulses  of  the  restituting  collision  (px  *,&>*)  * 


(Pl»  Pa). 

The  collision  cross  section  I(g,©)  is  usually  defined  by: 


-  Kie)Jn 

using  the  differential  solid  angle  dfl  and  the  scattering  angle  ©.  There¬ 
fore  one  obtains: 

-fVMA)  -pfc  jJaf  (3.35) 


Cto  substituting  Ofi»(Px  Inx^a)  into  (3. 31)  (in  place  of  Q*(xx  |^i^)],  and 
integrating  over  and  r^: 


(3-36) 


which  Is  exactly  the  collision  integral  in  the  Boltsaann  equation. 
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Turning  now  to  the  A2(x1|F1),  one  can  write  in  the  spatially  uni¬ 
form  case; 


(3-57) 


where 

-$>*>  |  kfl-p  (3-38) 

and  all  the  operators  are  omitted  since  they  have  no  effect.  Call 

the  x3-integral  in  (3.36)  0(x1X2),  then  one  sees  by  the  reasoning  used 

in  Section  4,  that  0(x1X2)  o  if  |r21  j  »  |  &2r0»  Furthermore 

it  is  not  difficult  to  show  by  the  use  of  centre  of  mass  and  relative  co¬ 
ordinates,  that  0(xx,xa)  depends  only  on  plf  p2  and  r21.  Therefore 
also  si^(xxx2)  0(xxx2)  will  depend  o^ly  on  these  variables,  and  for 
r  -*■  oo  the  result  will  be  zero  since  the  operator  will  completely 

separate  the  particles  1  and  2.  Next,  according  to  the  definition  of 
S  ^'(Xxxa): 


Calling  for  a  aoaent  S^(xxxx)  0(xixa)  «  ?^(x1x2 ) ,  then  one  can  write 
(5.38)  in  the  forms 

A  “  &k<  P&&*’ +  pt  &(w,v> 
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In  the  third  term  the  T-integral  can  be  done  and  since  the  upper  limit 
does  not  contribute,  one  obtains: 

0(*,v 

9  el*)  [ ( fts+AjSV’W-S™™ i (3.39) 

*1 

which  puts  in  evidence  that  Ha  depends  only  on  pj. .  Hence  also  Aa 
will  depend  only  on  pl#  as  it  should  in  the  spatially  uniform  case. 

6.  The  equilibrium  state 

One  should  expect  that  the  basic  kinetic  equation  (B)  will  lead  to 
the  state  of  thermodynamic  equilibrium  and  that  then  Fx  and  all  the 
higher  distribution  functions  will  agree  with  the  results  obtained  from 
the  mlcrocanonlcal  ensemble,  which  for  a  large  system  is  equivalent  to 
the  canonical  ensemble: 

£>"*(%••■  r*>  -A  *ff{-  -  feji-'M  1  (j.to) 


where  0  «  kT  and  A  a  normalization  factor.  From  (J.kO)  the  contracted 
distribution  functions  are  formed  as  before,  and  one  obtains  especially  for 


Fi 


the  Maxwell  distribution: 


(3.M) 


We  will  show  that  by  the  substitution  of  ple^(xi)  for  Fi(xiJt)  in  the 
kinetic  equation  all  the  functionals  A,(xx|F*)  became  identically  zero, 
and  that  for  the  higher  distribution  functions  one  obtains  virial  expan* 

slons  which  are  in  agreement  with  the  results  obtained  by  de  Boer1®  and 

19 


by  Mayer  and  Nontroll 
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a.  Zeroth  order. — Since 
from  (3.11): 


is  spatially  uniform  one  obtains 


inti? sc^5 


i«»  * 


Since  S are  the  momenta  of  the  s  particles  before  the 
s -tuple  collision  the  Hamiltonian  for  the  s  particles  can  be  written  as: 


Hence  it  follows  from  (3*^2)  that: 

ro)  1 


(5.M) 


an  to  be  expected.  From  (3.3t)  one  then  gets: 

Introducing  the  relative  coordinate  r  *  the  space  part  of  the 

x2-integral  becomes: 


which  is  clearly  aero,  Bence  Ai(xx|Fx)  -  0,  as  follows  also  immediately 
from  the  Boltzmann  form  (3*56)  for  Ax. 

b.  First  order.— Since  Ax  *  0,  Dxf|0^  ■  0,  and  since  Fx  is  spa¬ 
tially  uniform,  one  obtains  from  (3*18)  and  (3*^5 )  immediately: 

foi'lfi")  -  Jr  ($4+6.,)  e* 
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From  ^Hg,  exp  (-Hg/©)J  *  0  one  obtains: 

ewe* 


Introducing  this  in  (3. 1*5)  the  last  tera  in  (3.1*6)  clearly 
zero.  Since  as  we  saw  in  Section  5# 


integrates  to 


(3.1*60 


the  t -integral  in  (3.1*5)  gives: 


[i -£>.*;) 


Writing  H3  *  Ha  +  ~  p§  +  0i3  +  023#  and.  integrating  over  p3,  one  gets: 

2nt 


F^-'FV  e"^  fj£  e 


(5.1*7) 


Following  Mayer,  we  introduce : 


(3.1*6) 


then: 

Because  fjj  ♦  0  if  rjj  >  r0,  the  last  three  teras  are  constants  while 
the  first  tern  is  a  function  of  |$a-4i  I  which  will  be  zero  if  (qg-qj  |  > 
2 r0  so  that: 


•  O 


Bence  one  obtains: 
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*  (xrnjf  Cffi  fX 

*  fit,  l„L,  (3.1*9)  . 

from  vhich  follows: 

Ax^O  -  / k  4 FT)  (3-50) 

Since  according  to  (3.^9)  depends  on  the  coordinates  only  through 

r  a  |  it  follows  again  by  introducing  relative  coordinates  that 

the  space  part  of  the  Xa-integral  in  (3*50)  vanishes.  Hence  Aa(xx  |Fx  ^ )  * 

0. 

c.  Second  order.— To  find  one  needs  F^^XxXaXalFi*^), 

for  vhich  one  finds  analogous  to  (3*^7): 

c  ^  *’$“***  (3-51) 

Since  Ax  a  A*  a  o,  a  DiF^  ^  a  o,  and  from  the  spatial  uniform¬ 

ity  of  F^e  one  obtains  from  (3*19)  immediately: 

(8,,+ft,)  £%,*.*,  I (3-») 

In  order  to  get  an  explicit  expression  it  is  more  convenient  to  derive  a 
differential  equation  for  t Operating  on  the  left  side  of  (3*52) 


with: 
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using  the  operator  identity  (3. 46a)  and  carrying  out  the  r- integral,  one 
obtains : 


(3.53) 


since  it  is  easy  to  show  that  the  tqpper  liait  t  -  «>  gives  no  contribu¬ 
tion.  Following  Glauberman,^0  we  try  to  solve  (3.53)  by  the  Ansatz: 


The  left  side  of  (3-53)  becomes: 


(3  54) 


while. 


using  (3*51)#  the  right  hand  side  can  be  written  in  the  form: 


-  “fit  e  ^  ] 


where  a  is  the  constant  defined  by: 


Hence  from  (3-53)  one  finds  Immediately  the  particular  solution: 
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MV- i  Ipi/U  c  -**&&) 

-*rk%  (LUU+tMUUAt* 

+ UU+UU+UU,*UC+  L.M  (J.55) 

From  this  fora  one  sees  that  XeC^.,^)  depends  only  on  r  »  Itg-Qi  I 
and  it  goes  to  zero  if  r  -*•  »„  Since  according  to  the  general  boundary 
condition  (C±),  S^)(XiXs!)42)  must  be  zero,  equation  (3*33)  is  the 

solution  of  (3.53)  which  is  required.  Introducing  (3*55)  ia  (3.5M  one 
has  the  explicit  fora  for  F^^XiXalpI*^),  which  is  in  agreement  with 
the  result  of  de  Boer,  and  Mayer  and  Montroll.  Since  f|2^  depends  on 
the  coordinates  only  through  r  ■  |q2-qj.|,  one  proves  as  before  that 
*3(*ll^e))  »  0. 

It  is  clear  that  in  this  way  one  can  go  cm.  One  will  obtain  in  any 
orders 

where  X^(qila)  will  depend  only  on  r  ■  |  and  can  be  expressed 

as  integrals  over  combinations  of  Mayer  functions  tj  .  Hence  also  in  any 
order  At(xx  |f!**)  -  0. 


CHAPTER  IV 


THE  MACROSCOPIC  EQUATIONS 
PREPARATION  FOR  THE  u -EXPANSION 

1.  The  macroscopic  quantities 

The  usual  macroscopic  quantities  describing  the  state  of  the  gas  are 
obtained  from  the  first  and  second  distribution  function  by  further  aver¬ 
aging  over  the  impulse  variables.  They  are  defined  as  follows: 

a.  Number  density  n($,t). —This  is  defined  as  the  average  number  of 
molecules  in  the  volume  element  dq,  or  as  the  product  of  the  total  number 
of  molecules  in  a  system  and  the  probability  of  finding  a  particular  mole¬ 
cule  in  dif.  Or: 

In  the  following  it  is  sometimes  convenient  to  introduce  a  dimensionless 
function  v(q,t)  by: 

-ppft  (k.2) 

For  a  spatially  uniform  system  clearly  v  »  1. 

b.  Itocroscoplc  flow  velocity  u(&,t).  —This  is  defined  by  stating 
that  mil  •  nd$  is  the  average  momentum  of  the  molecules  in  the  volume  d$. 
Therefore: 


Jk 
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so  that: 

(4.5) 

c.  The  kinetic  temperature  6($,t ) .<— This  Is  defined  In  terms  of  the 
average  kinetic  energy  of  the  molecules  in  dq,  which  is  given  by: 


This  can  be  split  into  two  parts: 

The  first  part  is  the  kinetic  energy  of  mass  motion  and  the  second  part 
la  the  energy  of  the  random  motion.  We  now  define  ©($,t)  by: 

c.4) 

Clearly  in  equilibrium  o/k  (k  ■  Boltzmann  constant)  will  become  the  abao- 
lute  thermodynamic  temperature. 

d.  Internal  energy  density  n($«t)  c($»t).— This  is  defined  as  the 
sum  cf  the  random  motion  part  of  the  kinetic  energy  ^  n($,t)  0($,t)  and 
of  the  average  intermole cular  potential  energy  n($,t)  c^($#t)  which  is 
given  by: 

Bence: 
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2.  Tte  general  Macroscopic  equations 

The  equations  which  the  Macroscopic  quantities  satisfy  are  obtained 
from  the  first  two  of  the  B-B-G-K-Y  equations  (2,10)  and  (2.11)  by  aver¬ 
aging  over  the  impulse  variables.  In  these  equations  we  will  write  (q,p) 
for  (qx>Pi)#  apd.  we  will  use  the  tensor  notation  with  the  usual  summation 
convention.  We  will  use  Greek  letters  for  dummy  indices. 

a.  The  equation  of  continuity. —Integrating  (2.10)  over  $  one  ob¬ 
tains 

+  (»•«) 


By  dividing  by  v,  this  becomes  the  familiar  equation  of  continuity: 

’foj'-  -o  (*>.7) 


Introducing  the  "substantial  tine  derivative" 

&  •  £  ♦ 

this  can  also  be  written  in  the  form: 


(4.8) 


b.  The  general  hydrodynaadcal  equation.  —Multiplying  (2.10)  by 


and  Integrating  over  cat  obtains: 


where 


m+lf  ■M-fi’df.*) 


(4.10) 
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where  we  have  put  for  the  momentum  of  the  thermal  motion: 


A  ’ 

Pij(q,t)  is  the  familiar  expression  of  the  pressure  tensor  due  to  the 
kinetic  motion  of  molecules.  Using  (4.6),  equation  (4.9)  can  be  written 


in  the  form: 


nm  §?  * -  ^  -  Wh*l 


3RT 


(4.11) 


Now  it  can  be  shown  that 


(4.12) 


where  P^j($,t)  is  defined  by: 


Pifil  t)  *  ^ %;i) {U,15) 

o  s 


Therefore  (4.11)  becomes: 


(4.14) 


where: 


(4.15) 


is  the  total  stress  tensor.  From  (4,10)  and  (4.15)  one  sees  that  P1j($,t) 
is  symmetric. 

Proof  of  equation  (4.12)  (see  also  Enskog21):  Define  the  pair  den¬ 
sity  distribution: 

£<*«.*) 


38 


then  n2(|#q2*t)  is  a  syanetrlc  function  of  the  two  points  q  and  <ta,  and: 

Jw  i,  -  pi  k  Mi/)  **<t,  to 

where  ki  ■  (qi-q2i)/|$-4e|*  Now: 

?+Jx) 

Since  ng  is  a  function  of  q  ♦  &,  and  tr, 

ninx  l  w*, 

-VK 


Consequently 


2Ss*  -  i  iZut‘  jjJi  k  fa)  ( -  ntf-U,  f )} 


ijt  ki 


) 


using  the  sysnetry  of  n* ($,$•)* 

£•  The  energy  transport  equations .  —Multiplying  (2.10)  by  , 

integrating  over  p  and  rearranging  terns,  one  obtains  the  transport 
equation  of  the  kinetic  energy: 

*■ ^ 4  «•+  $ -  5-Ml.lt &*«.*> 

where  the  kinetic  part  of  the  heat  current  density  Is  given  by: 

‘  vf<l?4-£rK<1.lU  <*•*> 

Using  (k.lk)  and  (k.15),  one  can  also  write: 


(‘••17) 
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where  the  rate  of  the  deformation  tensor  Djj  is  defined  by: 

Di;  ('•.18) 

In  order  to  obtain  the  transport  equation  of  the  potential  energy 
density,  we  have  to  use  the  second  B-B-G-K-Y  equation  (2.11).  By  multi¬ 
plying  with  0( j$-4g| )/2vz,  integrating  over  p,  p2,  q-2  and  rearranging 
terms,  one  obtains: 

where: 

J-  (k-20> 

The  transport  equation  of  the  internal  energy  density  is  obtained 
by  adding  (4.17)  and  (4.19).  Ibis  gives: 

n#+i<iK+J& 

'irjtfWl  $£('&  +  '&)ECK*"'t)  ('••2D 

Similarly  to  (4.12),  it  can  be  shown  that  if  one  defines: 

J;  ’  £j[2u#4 (k-22) 


then 


(4.23) 
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and  considering  (4.12): 

where 

JififA)  -  & M*W  • 
jm+Jn1)  Ft  P, 


Thus  (4.21)  becomes : 


71§F  +  1£ 

where 


(4.25) 


(4.26) 


Note  that  the  potential  part  of  the  heat  current  density  (given  by  Jf  ♦ 
J^2)  depends  only  on  the  thermal  momenta. 

The  pressure  tensor  P^j(^t)  and  the  heat  current  vector  J^(<l»t) 

12d 

are  identical  with  those  of  Irving  and  Kirkwood. 

1*  Expansions  of  f£°  WxgjFx)  «£d  Ai(xi,  iFi) 

According  to  the  expression  (5.24)  for  Ax(xi|Fi)  and  the  formulae 
derived  in  the  preceding  section,  one  needs  only  to  consider  the  second 
distribution  function  for  phases  xx,  Xs  »uch  that  the  spatial  distance 
|5i-S*i  £  rQ.  We  shall  see  that  in  the  next  stage  (the  hydrodynamical 
stage)  of  the  approach  to  equilibria  all  distribution  functions  will  de¬ 
pend  oo  the  coordinates  through  sacroscopic  quantities  which  vary  slowly 


4l 

wjLth  position.  As  a  preparation  we  will  therefore  in  this  chapter  already 
expand  the  second  distribution  functions  in  powers  of  (r0  *  grad|  ),  and 
we  will  retain  terms  up  to  the  second  order. 

Since  according  to  (3*23): 

fihzjfr  ’{fa A  £«,  t)£«.  t>  Kx«  **  i  <  v  ik-£7) 

where: 

“  faZ-W-  Wl-  Oj'fsZSZ’Sct-IM-  Vi  (1,-28) 

we  first  expand  the  kernel  Kg.  Using  the  center  of  mass  and  relative  co¬ 
ordinates  : 

VI  -ft.  .  t- !•*««-** 

and  the  corresponding  Momenta: 

^  “R  *  ft~f*  m  9  ~  *fn  (*»*29b) 

one  shows  easily  that: 

ft.  “  <l,-50l>) 


so  that  the  operator: 
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S?T wSr'ku  -  e*p «p  (r  ^-J^) 

*  J^r  (^.31) 

depends  coaly  on  the  relative  coordinate  and  impulse  and  is  invariant  un¬ 
der  a  combined  orthogonal  transformation  of  r12  and  £12. 

Call: 

MAM-V  -7T.fM.EiH)  <*•*> 

putting  in  evidence  that  it  depends  on  the  spatial  coordinates  only  through 
r12.  One  further  has: 

sis -  Jj(Vbtm-T-{) 

(U.J5) 

Since  we  are  interested  in  those  values  of  r12  which  are  of  order  rQ  and 
since  for  such  rxa  one  readily  shows  that  \*im  ri2 1  is  also  of  order 
r0,  we  expend  (4.33)  in  a  Iky lor  series  around  .  One  thus  obtains: 


(frijTT.  + 


where: 


Tt'.dA.K  U)  =  - 

7TXi .  /h,  ?'.R  1 1.5!.)  “  7T, 


(M5a) 

(^•55b) 

(^•55c) 


*(^.35d) 


(^.35e) 


All  the  IT's  are  functions  of  ?i2,  ?2-^2  and  $ist,  which  have 

the  following  properties,  as  can  easily  be  verified: 

a)  Adding  a  constant  vector  to  all  momenta  will  not  change  the  TT's 
(Galilei  invariance ) . 

b)  Under  an  arbitrary  (proper  or  improper)  orthogonal  transformation 
applied  both  to  r12  and  all  the  momenta,  the  TP  s  transform  like  isotropic 
tensors. 

2)  For  an  interchange  of  the  two  particles  (corresponding  to  the 
transformation  (ria,  px,  p*,  4#  Hz)  ♦  (-?xa>  Pa#  Pi#  5a#  5i),TT^  and 
TTjt's  remain  unchanged,  while  theTT^'s  change  sign. 

Introducing  the  expansion  (k.Jk)  in  (k.2f)  and  carrying  out  the  in¬ 
tegrals  over  ti  and  leads  to  the  expansion: 


r<^iF)  -  c+ 


(*►-36) 
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where 


<£\  id  -  1  ^*} 


(4.37a) 


(4 .3Tb) 


(4.37c) 


All  the  j'  ■  are  to  be  considered  as  functions  of  <b,  rxa  and  p1#  pa. 
Because  of  property  (c)  of  the  functions  TT,  then's  with  the  upper  index 
A  are  antisyawtrlc  voider  an  Interchange  of  the  two  particles,  while  those 
with  upper  index  S  are  syanetric.  Bote  also,  that  because  of  (4.33)  the 
ant  isyaeetr  lc ^ •  can  be  written  as  divergences: 

C  ■il-5'C) 

C*  •  tt't'lT} 

Finally,  according  to  (3*24): 

/Vr-iP>  -  v 


(u.ja.) 

(V.») 
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with: 

jQ.rt  I  <,y  *  Jjr,  K,( *■  *  I  <.  <,)  C*.ko) 

Substituting  the  expansion  (4.34)  for  the  kernel  K2,  one  obtains  an  ex¬ 
pansion  for  the  binary  collision  operator: 

A,  -  a?  *  a"’ A+  <?;"*+  a“A + ars+ .  ^au 


We  omit  the  explicit  expressions,  since  they  are  quite  similar  to  the  ex¬ 
pressions  (4.37)  for  the  ^'s,  and  in  fact  are  obtained  from  these  by  op¬ 
erating  with  0x2  and  integrating  over  p2  and  rx2.  We  only  note  that: 


(4.42a) 


with: 


(4.42b) 


This  has  precisely  the  sane  fora  as  in  the  spatial  uniform  case,  which  was 
discussed  in  Chapter  III,  Section  5.  One  can  therefore  transform  #x°^S 
to  the  familiar  Boltzmann  form: 


k.  of  rP>ton»  Ift.)  jnd  *.txilFt) 

In  the  same  way  as  in  the  previous  section,  we  will  expand  F*(xxx2 |Fx ) 
and  A*(*i  )  *  These  functionals  are  connected  with  the  effects  of  triple 
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collisions,  and  as  a  result  the  formulae  are  rather  complex.  We  will 
therefore  only  indicate  the  method,  and  show  that  the  essential  proper¬ 
ties  of  the  previous  expansions  appear  here  again. 

According  to  (3-27),  one  can  writes 

m  Jjf.Js.kjf K/WlWJ  (4.W0 

with 

-  v  (M5) 

We  begin  again  with  the  expansion  of  the  kernel.  We  are  interested  only 
in  such  phases  xx,  x2  so  that  Iqi-qal  £  rQ.  In  Chapter  3,  Section  4, 
we  saw  that  for  such  phases  the  T-integral  extends  only  over  a  time  of 
the  order  of  the  collision  time  tq.  Also  in  the  integrand  of  Ka  only  such 
values  of  the  phase  x3  play  a  role  for  which  |qj.-qL3  I  and  | <12-^3  |  are 
of  order  r0.  It  is  therefore  convenient  to  introduce  center  of  mass  and 
relative  coordinates  according  to: 

*  £  (f.+l  :  VH  i  *  fr 1 

and  the  corresponding  momenta: 

R  -VP.+ F,  !  -  l-w-v 


(4.46a) 


(4.46b) 


Then: 


1*7 


and  one  finds : 

ss>as?«  ■  » *  l-rhkk+i&k~  1 '  4 '  &  fl 

S  ^  (‘•.‘•8a) 


(4.48b) 


(4. 1*  8c) 


Clearly  the  opeiators  J  depend  only  on  the  relative  quantities. 

Splitting  in  the  kernel  the  product  of  the  8-functions  in  the  coor¬ 

dinate  and  impulse  part  and  using  (4.48),  one  can  write  for  the  part  of 


the  integrand  in  K3  which  comes  after  the  ©-operators  (for  k 

JZsft 

•  { |  h%*  «£  S(ji- 1)} 


1,2,5): 


(*»A9) 
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where  R2  is  again  ^  (^l+qa)  and: 

*>»  *  I'  -  (f*>- &,)*£<&-&+ £-£H,-V 
X>»  -I'-  (|  V* 

X>*»-  (M.)  - 


One  can  easily  see  that  for  all  phases  where  the  three  distances  Iqi-q^l, 

Iqi-^l  and  {q2-q3 |  are  of  order  rc,  also  all  the  are  of  order 

rQ.  One  also  easily  verifies  that  by  an  interchange  of  the  particles  1 
Ik) 

and  2  the  change  according  to 


(J.f.S) 


*  %**(/> >.3)  ^  ’ 51 ) 

*.•  ))  *'V**(l.9.3'i 

^  '  '3 


One  now  can  start  the  expansion.  Just  as  in  the  previous  section, 
we  leave  the  impulse  part  ^  ^  1T&(Pj  -V  in  (4.49)  (which  depends  only 
on  relative  quantities)  as  it  is,  and  expand  the  coordinate  part.  Since 

*•  /  V  ) 

'  are  of  order  rQ,  we  write: 


(4.52) 


where  we  use  the  suamation  convention  both  for  the  Latin  index  1  as  for 
the  (Greek)  index  a,  which  denotes  the  vector  components.  Introduce 


/  k ) 

this  expansion  [extended  to  second  order  in  the  X^  ]  in  the  kernel  K3. 
Since  the  product  TT5($2+?3  r/j m-^)  commutes  with  the  ©£,  there  will 
occur  exparsions  in  which  the  remaining  S -operators  in  K3  act  on  thi prod 
uct.  It  is  easily  seen  that  one  can  write; 


=  J  z(i-^ \) 


(^.53) 


where 

pfl,»,3)  «  (?l»  £  r  (‘*•5*') 

Since  r  is  of  order  t0,  v  is  of  order  rc,  and  one  can  therefore  fi¬ 
nally  expand  in  powers  of  (-l/2  r12  +  v)  around  q^.  Doing  this  also 
up  to  second  order,  one  finally  obtains  the  expansion: 

K3(*,x.!<,5,y  >  JT Z<fr{)V0 

*  V’M 

+  "Hu {  % 


k 

where  we  have  introduced  for  abbreviation  the  differential  operator  7 
which  when  acting  on  a  product  of  functions  of  is  defined  bys 
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Thr  form  of  (4.55)  is  completely  analogous  to  the  expansion  (4.54)  for  Kg 
and  the  functions  cn  are  analogous  to  the  functions  TT.  They  are  defined 
by: 

‘  l Jh <ttj.{(b)+  9ai) X  ~ 

-Jue.X+e.a'l))*  scpt  -  %) 

=  fef»>  \k  I  (wS!"-‘>  -  s>>  fagfM + 

yo  y  1  (1*  55a) 

CO.u(\,}i.ft I’l.li.lp  2.  ^0  , 


(4.56b) 


0.Xy»£')J + 

+ «?!!'+ - 1)  (k.. 


(4,56d) 
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+  fti.»  -X  (0„  £*  &,  X’)] 


iia».  ^[(d^djA^'^Jt!-^  (8X'P>)X“+ 

+  (f  Wy)l(0V+BJJ?-jL(B.tJ?*  0.X>]}  . 

•  ft  (l‘-56e) 

In  here  the  only  undefined  symbol  is  the  operator  o^)  which  is  part  of 
the  two  particle  operator  S^(l,2)  expressed  in  the  two  particle  center 
of  mass  and  relative  variables.  One  has: 

St  <'•*>  -  **f  i  -r^Sr  ^  ]  •  °z'r 


It  is  easy  to  prove  that  the  functions  o>  have  the  same  properties  a), 
b)  and  c)  which  were  mentioned  on  page  1*3  for  the  functions  TT.  Intro¬ 
ducing  therefore  the  expansion  of  Ka  in  (4.44)  we  get  an  expansion  for 
P^(*i*a|*i)  which  is  of  the  same  form  as  the  expansion  for  Fg°^.  We 
write,  analogous  to  (4.36): 


•  •  e  •  * 


(4.58) 
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where: 


C<  'p>  4\W iftvivwt.uiu.v  (M9a) 

£/  *F7)  -{Ilmllfoi’V)  M-.UM.V  (U.59b) 

I  y/*!*1!})  (4.59c) 

%k/,pmw  «..**> 

£T<  if)  -  JfH^  \?{  pti .  V)  f  wv<» .».  > 

The  superscripts  S  end  A  refer  again  to  symmetry  or  antisymmetry  with  re¬ 
gard  to  an  interchange  of  the  particles  1  and  2.  One  still  has,  regarding 


and  ria  as  independent  variables  that: 


if-frCj 


(4.60a) 


r:  -  ki-*xr) 


(4.60b) 


Finally,  since  the  triple  collision  functional  is  given  by: 


A,»iF)  ft  v  (4.6l) 


with: 


XljfclW,)  - 


(**.62) 
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the  expansion  (4.53)  for  the  kernel  K3  leads  to  an  expansion: 


4<r,l  f->  -  Q[n‘  +  tff  +  +  Q™+Q?s- 


(4.63) 


similar  to  the  expansion  for  .  We  omit  again  the  explicit  expressions 

for  the  CX's  since  they  are  obtained  immediately  from  the  expansion  (4. 58) 
(1) 

of  the  Fa  . 


2*  Expansion  of  the  macroscopic  equations 

In  Section  2  we  have  derived  the  general  macroscopic  equations  from 
the  B-B-G-K-Y  hierarchy  of  equations.  These  equations  are  exact,  but  form 
only  a  general  scheme,  in  which  any  closed  (and  approximate)  system  of  equa¬ 
tions  for  the  five  macroscopic  quantities  will  have  to  fit.  In  this  sec¬ 
tion  we  will  derive  more  specific  macroscopic  equations  from  the  general 
kinetic  equation  (B)  of  Chapter  3,  in  which  we  will  use  already  the  ex¬ 
panded  forms  (U .Ul )  and  (4.63)  for  the  binary  and  ternary  collision  oper¬ 
ators.  The  basic  equation  Is  therefore: 


JF{ar<viE)+£(Q,;,*+Q:,i)}+ 


+  + .  <“•<*> 

Notice  that  this  equation,  Just  as  the  Boltzmann  equation,  is  a  differ- 
ential  equation  in  the  coordinate  q,  and  an  Integral  equation  in  the  mo¬ 
mentum  p , 

The  procedure  to  obtain  macroscopic  equations  is  the  same  as  used  in 
Section  2.  Integrating  (4.64)  over  p  and  dividing  by  v  gives  again  the 
continuity  equation: 


~w  +  &*«■)  -« 


since  the  right  hand  side  of  (4.64 j  does  not  contribute.  Then  by  multi¬ 
plying  with  p  and  integrating,  one  obtains  the  equations  of  motion  in  the 


form: 


7?7rt^^ 

Pt 


+ £  SJf  pi  d?+afr+~) 


+ £Pp  ^  {ar +ct‘*+  or*  ctM«y] 


(fc -65) 


Iv 

with  the  same  P^j  as  in  Section  2.  From  the  expression  of  the  (£'b 
in  terms  of  the  corresponding  expansion  (4.36)  and  (4.58)  of  the  second 
distribution  function  one  finds: 

Since  d^/br^  is  odd  in  r,  clearly  the  contribution  of  the  symmetric 
J^b  vanish.  Since  according  to  (4.38a,b)  and  (4.60a,b)  the  antisymmet¬ 
ric  F:  can  be  written  as  a  divergence,  one  getj  from  (4.65)  and  (4.66), 


the  equations  of  motion  in  the  desired  form: 


n  m 


J UL 


an 


(‘••67) 


with: 


Pj  *  P  f 

v  i# 


(4.68) 


mad) 


55 


Finally  by  multiplying  (4.61)  by  p2/^  and  integrating  one  gets  the 
kinetic  energy  equation  in  the  form: 

*  P*)  +  ^  +  + 

+  v*/^«  (4.70) 

Since ,  as  mentioned  at  the  end  of  Section  3>  (2i°^S(q#P  |FX )  can  be  trans¬ 
formed  into  the  familiar  Boltzmann  form,  its  contribution  in  (4.70)  will 
be  zero  according  to  a  familiar  argument  (see  Chapman  and  Cowling,  p.  67). 
In  general  one  has: 


and  therefore  in  general  both  the  symmetric  and  antisymmetric  ^ '  s  will 
contribute.  Splitting  p^  according  to: 

f>.  *  ft  ~  fa-  +  ft xr  ^  (4.72) 


then  the  second  pert  which  is  ayaaetrlc  In  the  two  particles  together 
with  the  antisymmetric  jF*§  can  be  written  as  a  divergence .  One  finds 
for  this  part: 


56 


with: 


The  first  part  of  ^4.72)  combines  with  the  symmetric  ¥  's,  and  altogether 
the  kinetic  energy  equation  becomes: 

«&(#*)  •  -  fcl'+Jt')  +  <*•*> 

with: 


_X 

V1 


'>C*C4 . 


Since  in  the  kinetic  stage  the  equation  for  F2(xxx2|Fx)  has  the 
same  fora  as  the  second  B-B-G-K-Y  equation,  the  equations  for  the  poten¬ 
tial  energy  will  be  the  same  c.s  (4.19)>  namely: 


(4.76) 


with: 

afityp  -  ifrllfW >F>  (V-TT) 


jftfiF)  -  M  I K ) 


(^.78) 
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Using  the  expanded  forms  (4. 56)  (U . 56 )  for  and  Fg1^,  clearly  the 

right  hand  side  of  (4.76)  becomes  -R(q|Fi).  In  only  the  symmetric 
jf*  '8  will  appear,  while  in  J*2  because  of  the  factor  jpj/m  both  the 
symmetric  and  the  antisymmetric  y<B  will  contribute.  Adding  (4.74) 
and  (4.76)  gives  the  total  energy  equation  in  the  same  form  as  in  Section 
2. 

What  is  achieved  is  the  expression  of  the  stress  tensor  PAj  and 
the  total  heat  flux  vector  Ji  in  terms  of  the  first  distribution  func¬ 
tion  only.  This  will  be  the  starting  form  of  the  macroscopic  equations 
in  the  next  hydrodynamic  stage. 


CHAPTER  V 


THE  THEORY  OF  THE  HYDRODYNAMIC AL  STAGE 

1.  The  basic  equations 

As  explained  in  the  introduction,  we  assume  that  as  the  gas  relaxes 
further  towards  the  equilibrium  state  a  second  coarse  graining  in  time 
(over  a  time  of  order  t0)  is  needed  in  order  to  describe  the  slow  varia¬ 
tion  of  the  macroscopic  quantities  in  time.  On  this  "hydrodynamic"  time 
scale,  we  assume  that  the  first  distribution  function  F*  depends  on  time 
only  through  the  macroscopic  quantities  n,  u  and  9,  so  that: 

R(*:t)  - ‘E'tf.fl  n>  (D) 

where  as  before  the  vertical  bar  denotes  that  F*  depends  functionally  on  n, 
and  0,  which  contain  the  whole  time  dependence.  The  form  (D )  is  assumed 
to  be  valid  for  any  initial  distribution  Fi(xjo)  after  an  initial  period  of 
order  t0.  The  macroscopic  quantities  are  expected  to  vary  smoothly  on  the 
hydrodynamic  time  scale,  and  to  fulfill  the  basic  hydrodynamic  equations  of 
the  form: 

(f  ) 

*  U  (f  |  n,  6) 

*  (!)  (f  I  U,  0) 

The  unknown  functionals  Fi($,  $Jn,  3,  9),  u.  9)tU($|n,  5,  9)  and 
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ettln,  «,  9)  must  follow  from  the  kinetic  equation  (^.64)  and  the  correopond- 
ing  macroscopic  equations  (Chapter  IV,.  Section  5)-  They  will  be  determined 
again  in  successive  approximation.  Since  at  this  stage  the  deviation  from 
equilibrium  is  caused  by  the  non-uniformity  of  the  macroscopic  quantities, 
the  development  parameter  (called  p)  will  be  a  measure  of  the  spatial  varia¬ 
tion  of  the  macroscopic  quantities.  Physically  u  will  be  of  the  order  of 
the  relative  variation  of  the  macroscopic  quantities  over  a  mean  free  path. 
However,  it  is  more  convenient  mathematically  to  use  p  as  a  formal  unifor¬ 
mity  parameter ,  with  u  »  0  corresponding  to  the  completely  uniform  state. 

The  expansion  of  the  basic  hydrodynamic  equations  (F)  will  .hen  be  of  the 
form: 

if  ■  ft . 

If-  *<r* 


/'*/”+ . 

(5-la) 

. 

(5-lb) 

. 

(5.1c) 

and  the  expansion  of  Fi  will  then  be  of  the  form: 

n,d,&)  + .  (5-2) 

Note  that  (5.1)  gives  the  time  derivatives  of  n,  u,  0;  these  quantities  them¬ 
selves  are  still  determined  from  Fi  in  the  usual  way  for  all  values  of  u.  We 
require  therefore  that: 


*<!.«>  -  P?  (. 

(5-5a) 

-  Pf  *  /. 

(5- 5b) 

(5-5c) 

6o 


and  for  i  *  1,2, 


jdf  ti 


!Jf  ir{i 


j  rf  li 


-  0 


(5.4) 


2.  Comparison  with  the  macroscopic  equations 

We  will  now  compare  the  hydrodynamic  equations  (5-1)  with  the  macro¬ 
scopic  equations  derived  in  Section  5  of  Chapter  IV,  in  which  only  the  first 
distribution  function  occurs.  Putting  in  the  expansion  (5*2)  will  lead  to 
a  n -expansion  of  the  macroscopic  equations,  if  in  addition  every  differentia¬ 
tion  after  the  coordinate  ia  multiplied  by  the  uniformity  parameter  u.  This 
is  necessary  for  consistence  since  the  differentiation  by  $  operates  only  on 
the  macroscopic  quantities  and  is  therefore  just  a  measure  of  the  non-uni - 
formitj  of  the  gas. 

To  abbreviate  the  formula  we  will  use  the  following  notation.  In  the 
macroscopic  equations  of  Chapter  IV,  the  functionals  J?  of  Fi  occur,  which 
involve  Fi  as  a  product  such  as  ,7 T  Fx($,tJi).  Therefore  introducing  (5*2), 
one  gets  in  zeroth  order  only  f0,  in  first  order  f0  and  tx,  in  second  order 
f0,  fi,  and  f 2,  etc.  We  write  therefore  the  ^-expansion  of  a  general  func¬ 
tional  Z(lFi)  of  this  type  in  the  form: 

ZdiK)  -  Zci L>  +  />Z (!(.,{,) 

+  /•{ztico  +  z  f  i 

With  this  notation,  one  obtains  f r  m  (4.68)  and  (4.69)  for  the  stress  ten¬ 
sor  the  expansion: 


(5.5) 


6i 


where : 


P*.i  -  P/t  +  P*i 


0,  I'  * 


(5.6) 


P*!t( 


lu°>  I'*’"'  (5- 7a) 


♦(OOETH*.' ••)) 


(5.7c) 


Analgously  one  obtains  from  (4.75)  and  (4.78)  for  the  heat  current  vector  the 
expansion: 

Ji  -  Ji,.  +/J..  v'J\*  +  •••  • 


(5-8) 


where : 


Ji.«  •  J.',  *t,  +J:*t  -  K” 

Ju<W>  -frlg-t, 


1*0.  i.  *.  •*  • 


■  -M?  llnyp*- 

■irM»z:<’h* . ) 


(5.11a) 
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'Itt'b'KaO* . ) 


(5.12b) 


# 

Finally  the  functional  R  defined  by  Eq.  (4.75)  can  be  expanded  in  the  form: 


R  «  ^  R,  ,  + - 


where : 


(5.15) 


(5-lM 


ra  i  M)  -  {%*( , C) ♦ 

.  +{f£'“U +  $?<&+-•}) 


(5.1*»b) 
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Equation  (5*13)  starts  with  the  first  power  of  u,  since: 

MtnXB)  * S#l(tiWsteilT£<ib+ . ) 

is  equal  to  zero,  if  fQ  is  the  Maxwell  distribution  in  p  (compare  also  Sec¬ 
tion  4). 

Substituting  in  the  macroscopic  equations  (4.67)  and  (4.74)  these  ex¬ 
pansions  for  Pjj*  and  R,  one  can  then  compare  the  continuity  equation  and 
these  equations  with  the  hydrodynamic  equations  (5.1).  Equating  equal  powers 
of  u  one  obtains: 


>f"V|  I' n.ZJ) 

^ nuu 

(3.15a) 

»  o  i  •  2.  3,  -  ' 

(5.15b) 

(jfffiw.a.fi) 

+  M  ^  J 

(5.16a) 

• 

• 

-  . JLRjLI 

nm 

(5.16b) 

• 

0ul{ln.3.$) 

(5.17a) 

(5.1Tb) 

• 

• 

* 
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The  internal  energy  equation  obtained  from  (4.74)  and  (U . 76 )  does  not 
play  any  particular  role  for  the  determination  of  the  hydrohynamie  equations, 
but  since  we  need  the  result  later,  we  will  also  expand  this  in  powers  of  p 
analogous  to  (5-l): 


(5.16) 


where : 


E°Vf  a>  -  -  +  ,„}-•*  'Iff 

(519b) 


and  €0,  €1...  are  the  result  of  the  p-expansion  for  the  internal  energy  e. 


2*  The  expansion  of  the  kinetic  equation 

So  far  the  development  is  purely  formal.  The  functionals  ^  and  8^  ^ 
in  the  hydrodynamic  equations  as  given  by  (5.16)  and  (5.17)  still  depend  on 
the  unknown  functionals  f0,  f ,,  etc.  These  ha/e  to  be  determined  from  the 
kinetic  equation  (4.64),  Just  as  in  the  kinetic  stage  the  unkrown  functionals 
Fa  had  to  be  determined  from  the  B-B-G-K-Y  equations. 

To  do  this,  we  first  expand  the  kinetic  equation  also  in  powers  of  p. 
Writing: 


+ 


/-rvr 


(5.20) 


I*  %  •*  ♦  •> 

where  the  p'*/all  depend  on  q,  p  and  functionally  on  n,  u  and  d,  one  obtains 
by  comparing  (5.20)  with  the  p-expansion  of  then's  (as  given  by  Eqs.  (4.4l) 
and  (4.63)  of  Chapter  IV,  using  (>.2)  and  again  multiplying  all  derivatives 
after  q  by  the  uniformity  parameter  p): 


65 


(5.21a) 


Hiring) 

+  W'H.IM’KO  <Jn{J,) + tf>U>+  -  • } 
+[CHJ>^aU.V---j 

(5.21c) 

.  +  {q?s< '  l  o + i  •  •  ■  ■  j 

Since  Fj.  depends  on  the  time  only  through  the  n,  $,  and  0,  one  haB  [similar 
to  (5«^a)]: 


where  the  6  denote  functional  derivatives  and  the  bn/bt ,  bu/bt ,  bQ/bt  are 
given  by  the  hydrodynamic  equations  (5-l)«  Introduce  now,  similar  to  the 
operators  in  the  kinetic  stage  (p.15),  the  operators  such  that  for 

an  arbitrary  functional  t(|n,uy9): 

-  [li-N'O+tli-O  ■*■[!}■  ^‘3  (5-22! 

Using  these  operators ,  expanding  the  lefthand  side  of  (5*20)  and  equating 
equal  powers  of  w,  one  obtains: 


66 


o  = 

(5.23a) 

sn  -  r'” ' 

(5.23b) 

vmi  +  5"’ t,  -  r* 

(5.23c) 

and  so  on. 

These  equations  are  the  basic  integral  equations  which  must  be  solved, 
subject  to  .the  auxiliary  conditions  (5*3a,b,c)  and  (5-M*  Note  that  the 
first  equation  (5 -23a)  involves  only  fQ  and  should  therefore  determine  f 
completely.  Knowing  f0  allows  one  to  find  the  hydrodynamic  equations  in 
first  approximation.  Therefore  the  second  equation  (5* 23b)  becomes  an  in- 
iegral  equatton  for  fx.  Knowing  fi  one  can  find  the  hydrodynamic  equations  in 
the  next  approximation,  and  in  this  way  the  successive  approximation  method 
goes  along.  We  will  only  discuss  the  first  two  approximations,  which  corres¬ 
pond  to  the  Euler  and  Stokes -Navier  forms  of  the  hydrodynamics 1  equations. 

4.  The  ideal  fluid  equations  (Euler ) 

Clearly  one  should  expect  that  the  solution  of  (3.23a)  for  f0  is  the 
local  equi librium  distribution: 

IU.9)  (5.2M 
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The  auxiliary  conditions  (5*3a,b,c)  are  then  clearly  satisfied.  It  will  be 
sufficient  for  our  purpose  to  show  that  (5.24)  satisfies  the  equation  J"^0^  0. 
From  the  properties  of  the  functions  «0  [see  (4.32)]  and  [see (4. 56a)  3  dis¬ 
cussed  on  p.  45  and  especially  from  their  Galilei  invariance,  one  easily  ob¬ 
tains  from  (4.37a)  and  (4.59a)  if  Fj.  is  given  by  (5.24),  that: 


®  »  o-  u. 


(5-25) 


These  pair  distribution  functions  have  therefore  the  same  form  as  in  equili¬ 
brium  except  that  the  f  and  ^  occur  instead  of  the  ^  and  p,  .  As  a  re¬ 
sult  one  proves  exactly  as  in  the  equlli-rium  case  (Chapter  IIT,  Sect.  6) 
that  all  thefl^°^S(q,p  |fQ)  are  zero,  so  that  also  *  0. 

There  remains  the  question  whether  (5.24)  is  the  only  solution  of  the 
equation  J"^°^  *  O.  Since,  as  we  saw,  dt°^B  can  be  transformed  to  the  fa¬ 
miliar  Boltzmann  collision  integral,  one  can  appeal  to  the  H-theorem  to 
show  that  (5.24)  is  the  only  function  which  makes  dj°^B  equal  to  zero. 

We  have  not  attempted  to  generalize  the  H-theorem  to  show  the  same  fact  for 
the  higher  ^'s. 

Using  the  explicit  form  (5*24)  for  f0,  from  (5«5)»  (5-6),  and  (5*7a,b) 
one  obtains  the  zeroth  approximation  of  the  pressure  tensor  in  the  form; 

1  *  r<? 


where  the  scalar  pressure  p  as  function  of  the  density  n  is  given  by  the 
equlllbriua  virial  expansion,  which  is: 

f  •  nfi  [  »  - ) 


(5.27) 
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where 

M)  •  ~h  {dr,  t, 

4«'  -tfaJlLU, 


using  as  in  Chapter  III  the  Mayer  notation: 


For  the  proof  of  (5*27)  see  Appendix  I. 

Using  (5*26)  one  gets  from  (5.16a): 

uyf  In, u, 6)  =  +  lml£) 


(5-28) 


and  therefore  in  this  approximation  the  hydrodynamic  equations  become  the 
Euler  equations: 

15.  ,  .  2f.7 LM 

dfc  (5.29a) 


7t 


(5.29b) 


Note  that  the  functional  e^($jn,u,©)  is  not  completely  determined  by 
fQ,  since  Ri  involves  f*.  In  this  approximation  it  is  therefore  not  possible 
to  write  down  an  equation  for  d©/dt.  However  it  ie  possible  to  find  an  energy 
equation.  First  note  that  from  (5.10),  (5.11a)  and  (5.12a)  follows  that: 


(5-50) 


As  to  be  expected,  the  total  heat  current  is  therefore  zero  in  this  approxi¬ 
mation.  From  (5* 19a)  one  then  finds  that: 
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e"’  ■  -(“■?£ 


(5.31) 


with: 


8.  =2 :6*£t 

~  i6  +  J^JPH  &  4  IK'”S(IU  +  K‘hi.)  +  --'  }  (5.32) 

--j-e  +  ijdU  C~£  {  n  +  n'ftl  fJo  + . } 

5  ^  f  i  . } 

The  energy  density  €q  is  therefore  found  from  the  virial  expansion  of  p 
according  to  the  formula  of  equilibrium  thermodynamics.  Hence  the  energy 


equation  becomes 

2t 


_  ,,  2£±. 

U*W.  **£ 


(5.33) 


expressing  the  fact  that  all  changes  occur  adiabatically  in  this  approximation. 

For  the  following  we  will  need  a  simplified  expression  for  8^,  which 
can  be  obtained  since  the  second  part  of  the  expansion  for  Ri  [see  Eq.  (5*l^a)J 
can  be  evaluated.  One  finds: 

where : 

c"  -  ii*M & t'Vl.V V* %,t &i ini)  »•*•) 
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(5.35b) 


Prccf:  By  explicit  calculation  one  finds: 


We  mentioned  in  Chapter  IV  that  the  it’s  and  u>'s  are  isotropic  tensor  fields 
depending  on  r,  M  ,  and  t  .  As  a  result  the  A >  B,  and  C  are  isotropic 
tensor  fields  depending  on  r, ^  and  $  .  For  a  general  isotropic  tensor 
*ik...l  of  this  *ind>  the  integral: 

vi 11  be  a  numerical  tensor  of  the  same  rank.*  Hence  the  numerical  tensors 


*For  a  proof  of  this  theorem,  and  some  information  about  isotropic  and  numeri¬ 
cal  tenrors  see  Appendix  II. 
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obtained  from  M*  .M  )  and  8(^jP,  f.)  must  be  zero,  while  the  numerical  ten¬ 
sor  obtained  from  the  C  ($,$,%)  must  be  constants  multiplied  by  the 
Kronecker  Therefore  the  second  part  of  R2  becomes: 


From  (8.1Tb)  using  (5*26)  and  (5-30)  one  then  finds  ( 5 • 3*0 • 


2‘  Detenalnav i on  of  the  form  of  f  2 

a.  General  method. --To  derive  the  hydrodynamic  equations  in  the  next 
approximation,  one  first  her.  to  determine  f2  from  the  kinetic  equation  (5.23b). 
From  the  definition  (5.2?)  of  the  operator  Dvl^,  one  has: 


=  L  {n  N"’+  Ui  +  (-^t  '  h-56) 

«  <n  Q'  +  *i*  Qt  + . 

•rfhere,  using  the  results  for  IJ and  0^*^  from  the  previous  section 

and  writing  again  fQ 


(5* 37a) 
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+  '  £-X-$-  -  fP)^ ft  (5.37b) 

-3rf  jfe-  £>/#  $  J tyJf&P-Tutot.  0] 

By  writing  *8 ^^fo,  except  for  so  far  as  fj.  is  concerned,  as  a  power  series 
in  the  density  n,  we  have  separated  the  density  dependence  due  to  the  inter¬ 
action  of  the  molecules  in  pairs,  triples,  etc.,  from  the  dependence  on  the 
relative  change  of  the  macroscopic  quantities  n,  u,  and  0,  which  are  a  meas¬ 
ure  of  the  non-uniformity  of  the  gas  and  which  are  taken  into  account  up  to 
the  first  order.  One  can  also  say,  that  the  uniformity  parameter  u  is  a 
function  of  n,  since  u  measures  the  relative  change  of  n,  IS,  and  0  over  a 
mean  free  path,  and  at  higher  densities  the  mean  free  path  X.  is  no  more  in¬ 
versely  proportional  to  n  and  should  be  written  in  the  form: 


(  /  Vo^T + 


(5.36) 


In^Eq.  (5.23b)  we  will  therefore  also  develop  the  functional  r^’ln 


powers  of  n.  One  obtains: 


P".-«  + 

^  ^  (5.39) 

+*•[<?;%.) + cer^v + oiuo 

^  . » . . . 

where  for  abbreviation  we  have  put: 

»t>-  twM  &«%•  (540) 


We  will  now  seek  a  solution  for  fj  also  in  the  fora  of  a  power  series 
in  n.  Putting:  • 

i,  -  Lc  +  ”  L  + .  VM) 

and  equating  in  (5«2Jb)  equal  powers  of  a,  one  obtains: 


q,  +*(%£•*+ -oratu 

WtU  -  Q^ov-arav  -Q'Z'atU 

-ctr\  i  to 


(5.42a) 


(5.42b) 


These  are  inhomogeneous  linear  integral  equations  for  f j.)0,  ?ifi  which  have 
to  be  solved  successively.  Since  the  auxiliary  conditions  (5.4)  must  hold 
for  all  densities,  both  fx,o  and  f  1,1  must  fulfill  these  conditions. 

b.  The  Chapman -Enskog  theory. --In  our  formulation  this  theory  gives  the 
solution  for  fj.,o*  The  left  hand  side  of  (5*42a)  can  be  simplified  and  written 
in  the  form: 


The  right  hand  side  is  a  linear  integral  operator,  which  is  iso¬ 

tropic  in  the  f -space.  One  now  can  aake  use  of  the  following  theorem: 

a  linear  isotropic  operator  in  ^ -space  and  Rij.,.»(^)  is  an 


isotropic  tensor  in  this  space,  then  a  solution  of  the  inhooogeneo'ir  equation: 


will  be  an  isctropic  tensor  of  the  seat  rank  as  R.  If  the  homogeneous  e- 
quation  J(f>  *  0  has  no  solutions  this  isotropic  tensor  will  be  the  only 
so] ut ion.  Otherwise  we  have  to  add  to  the  isotropic  tenser  a  linear  com¬ 
bination  of  the  solutions  of  the  homogeneous  equation  multiplied  with  the 
appropriate  numerical  tensors. 

For  an  indication  of  the  proof,  see  Appendix  II.  In  our  case  the 
homogeneous  equation ■  0  has  the  five  solutions  1,  and  corres¬ 
ponding  to  the  number,  momentum  and  energy  conservation  in  the  binary  col¬ 
lision.  Since  the  solution  of  (5-^2a)  has  to  be  orthogonal  to  these  five 
quantities  (with  weight  yo)  according  to  the  auxiliary  conditions  (5*1*),  it 
is  clear  that  the  solution  will  be  uniquely  determined.  Define  the  functions 
Vll0)(?)  and  «ij(o)(?)  as  the  solutions  of  the  following  integral  equations, 
which  are  orthogonal  to  1,  fa  and  with  weight  function yQ: 


Qrbf..%vn  •-$(£--£)£ 


(5.^a) 


&  Cl?.. Jwip  = 


According  to  the  theorem  V*  ^(p)  and  must  be  an  isotropic  vector 

reap,  tensor  field,  since  the  right  hand  side  of  (5*^3a,b)  have  this  property. 
Therefore  Vj^0^  and  wi^°^  must  have  the  form: 


Vfc?) 


V"(f)  -it- 


(5.^a) 


Wy  C?)  *  WlpC&f-  - 


(5.U4b) 


(o)  (o) 

where  V  and  V  are  two  scalar  functions  of  |f  j. 
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From  the  form  (5-^b)  one  sees  easily  that  Wj^0)  is  orthogonal  to  1,  and 
jp:  this  is  also  the  case  for  V^^°^if: 

j  df  -  0  (5.1*5) 

In  terms  of  these  functions  one  can  then  express  f\,o  in  the  form: 

I n‘tJ)  “  "  %(f) {  v'^A-W? 

u  **  (5-^6) 

f (M  -  } 

There  remains  the  problem  of  finding  the  scalar  functions  and  V^°K 

This  can  only  be  done  in  successive  approximation,  either  by  using  varia¬ 
tional  methods  or  by  developing  and  in  an  appropriate  set  of  ortho¬ 
gonal  functions.  For  the  details  see  the  book  of  Chapman  and  Cowling,  Chap¬ 
ter  8. 

c.  *ff»ct  of  triple  collisions.  Knowing  fx,o>  the  left  hand  side 

of  (5.^2b)  is  completely  known,  and  we  can  therefore  determine  fi,x  in  the 
same  way  as  fi,o*  Note  that  the  Integral  operator  in  (5»^2b)  is  the  same 
as  in  (5- 42a)  and  involves  therefore  only  the  binary  collision  cross  section. 
The  triple  collisions  enter  in  the  left  side  of  (5.42b)  because  of  the  term 

|Jo»fi,o)«  The  left  hand  side  therefore  becomes  much  more  complicated. 
Using  the  isotropic  tensor  property  of  and  the  properties  of  numeri¬ 
cal  tensors  one  can  show  that: 


The  proof  is  completely  analogous  to  the  calculation  of  the  second  part  of 
Ri  outlined  on  p.70.  That  the  result  is  zero  is  due  to  the  fact  that  since 
the  trace  of  Wjj ia  zero  the  corresponding  numerical  tensor  of  the  second 
rank  must  also  have  zero  trace  and  is  therefore  identically  zero.  Therefore 

(fat  becomes: 


* 
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Qx(u.,ti  -  y.  c  *  ft  ^ 

+6rly-4-Xit-iF)? frJ 

By  explicit  calculation  one  finds: 

*  7.  [«4>^  +  4<H£  + 

*  a  c  wNf* + W%j 

CCatU  ■  5.  C  + A/Hl 

where 

a,<?>  *  Pi 

“  l(j*j~  ^')|L 

Cq(f>  *  “  Sy) 


^V^"CFrW?* 


(5.^7) 

(5.W) 

(5^9a) 

(5^9b) 

(5.^9c) 

b 

(5.^9d) 

(5- ^9«) 

(5^9f) 

(5-*9*) 
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In  here  the  %h\t  and  co#  [given  by  (U.56a)  of  Chapter  IV)  are  writ  en 
as  functions  of  and  r.  Because  of  their  properties  mentioned  on  p.l*3i 
■§l,  "ij»  H  end  are  isotropic  tensors.  Furthermore  since  the  trace  of 


zero,  the  trace  of  jty  is  also  zero.  Therefore  one 

• 

can  write: 

«£<?>  “ 

(5.50a) 

Cij0)  =  C/f)(-4f--  -Jmg-fy)  +C,(f)£j 

(5.50b) 

(5.50c) 

Mi/f)  -  /*<(>(*$’  " 

(5-50d) 

Therefore  (5«**2b)  can  be  written  in  the  form: 

-  %  [UV^  +  %<$%}  (5-51) 

with: 

Li<b  *  [ i<P‘f  -(-£f  -  {■)$■  - 

(5.52a) 

M«Cf)  « [-(f  +/<r>xty  -  3S  fy) 

(5.52b) 

+ BFCi-r  *  i)  +  **<v}  Stj  ] 

Rote  that  the  ter*  with  d  log  n/dq*  cancels. 

Corresponding  to  the  V^°)(^)  and  we  introduce  V^1^)  and 

as  the  solution  of  the  integral  equations: 
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arc  it.  jo  » 


(5.55a) 


arc  1 


(5.55b) 


One  can  conclude  again  that  the  and  Wj.j^  are  isotropic  tensor  fields 


and  that  they  must  have  therefore  the  fora: 


V?'(t)  -  Y"(f)  ■*£ 


=  W7’ff  )(-£$-  -  jlj  Sij)  +  W^pSij 


(5.5^a) 


(5- 5^b) 


where  V^1^,  Wx^  and  W2^  are  three  scalar  functions,  which  still  must  ful¬ 


fill  the  conditions: 


(5-55) 


in  order  to  fulfill  the  auxiliary  conditions  (5.*0-  One  gets  two  functions 
because  the  trace  of  W^j  ^  ^will  not  be  zero  since  the  trace  of  MiJ  is 
not  zero. 


6.  The  Stokes -Wavier  equations 

With  the  known  fora  of  the  distribution  function  f*  one  can  find  the 
stress  tensors  and  Pijfi  given  by  (5- 7a),  (5* 7c)  in  first  approxima- 

tion  and  hence  the  hydrodynamic  equations  up  to  order  u*.  The  stress  becosies: 


Ri»»  *  f g,»  +  ?l  * 


(556) 


where: 
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n,  -  r  +  •••■• 

v  «  + . 

(5.57) 

and: 

1 

C  =  %<r>  w‘*Vr> 

(5.58a) 

7 

I 

f,"  -IsItfT&f.wrf'w- 

-  irfaW,  $  vvvy  • 

•  t  w,oVi,){^- -Jg^) x<v. 

•  (5.58b) 

1 

*  m  £  3,(t>  h*  ( r> 

(5.58c) 

The  total  stress  tensor  p6ij  ♦pij,!  has  therefore  the  familiar 

Stokes- 

Navler  form  for  viscous  fluids.  Two  viscosity  coefficients  appear,  the  shear 
viscosity  coefficient  qx  and  the  bulk  viscosity  coefficient  for  which 
one  obtain*;  expansions  similar  to  the  virial  expansion  for  the  pressure.  The 
zeroth  approximation  is  the  Chapman -Enskog  value;  consists  of  two 

parts,  one  depending  on  the  triple  collisions  (which  come  in  through  Wx^) 
and  the  other  depending  on  the  potential  energy  0(r).  The  bulk  viscosity 
t|s  is  in  this  approximation  proportional  to  the  density  and  is  a  consequence 
of  higher  order  collisions  through 

The  hydrodynamic  equations  up  to  order  n*  become: 

$  •  - 

v  -  -(**♦  *  A  ^Jf ) 


(5-59) 


8o 


where  to  be  consistent  we  should  use  for  the  pressure  p  the  virial  expan¬ 
sion  up  to  the  third  virial  coefficient. 

By  substituting  the  fora  of  fi  into  the  formula  (5*10),  (5. lib),  (5.12b) 
feu  the  heat  current  density  and  taking  terms  up  to  order  fl,,  one  obtains: 

Ji  =  Jo  +  J t  +Ji?  (5-60) 

with: 

•X  =  Tf0)  +  n  T<n  + .  (5.61) 


and: 

T<*’  =  &-T/PW?) 


(5.62a) 


r<”  =  3~fJf  Tm*  %(P>  Vty 

-  %<vvv  • 

•  7£o?,p.  #  1  f,.?,)  (5-62b> 

+  4  %<vi<v  • 

•  ,  rc  v  ^"7?;  c*<  p  •  ?,  1  i  ?. ) 

Equation  (5.60)  is  the  Fourier  law  for  heat  conduction  and  r  (except  for  a 
factor  (k/9))is  the  heat  conductivity  coefficient,  for  which  one  obtains 
again  a  virial  like  expansion.  The  zeroth  approximation  t^° ^  is  again  the 
Chapman-Enskog  value.  The  first  approximation  consist  like  of  a 

triple  collision  and  a  potential  energy  part. 

As  In  Section  U,  it  is  not  possible  to  write  down  the  9-equation  up  to 
the  corresponding  order  in  (i,  since  0^  depends  on  fa  and  Is  therefore  not 
completely  determined.  From  (5*3*0  one  can  calculate  and  one  can 
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simplify  the  expression  for  8^*^.  However,  since  the  results  would  only  be 
needed  in  the  next  approximation,  we  omit  the  detailed  formula. 

One  can  write  down  the  energy  equation  up  to  order  u2.  The  result  is: 

+ 1 1  E,Af  +  -tff)  ^ 

with: 

l  =  £„  +  <sf 


where  <J#  is  given  by  (5*32)  and: 


* 


(5.65) 


CHAPTER  VI 


COMPARISON  WITH  THE  ENSKOG  THEORY  OF  DENSE  GASES 

1.  Introduction 

Q 

For  the  special  molecular  model  of  rigid  elastic  spheres,  Enskog0  de¬ 
veloped  many  years  ago  a  theory  of  dense  gases,  in  which  only  binary  col¬ 
lisions  were  considered.  For  an  account  of  this  t,.  .»ory  see  also  Chapter  16 
of  the  book  of  Chapman  and  Cowling.*^  In  order  to  compare  Enskog' s  theory 
with  ours,  we  will  calculate  the  first  density  corrections  to  the  viscosity 
and  heat  conduction  coefficients  from  our  formula  using  the  elastic  sphere 
model.  We  will  see,  that  the  results  are  identical  with  Enskog' s  results 

if  we  consider  only  those  parts  of  the  density  corrections  which  are  due  to 

"  ♦ 

the  binary  collisions.  However  we  were’ unable  to  calculate  explicitly  the 
contribution  of  the  ternary  collisions,  and  it  is  quite  unlikely  that  this 
contribution  vanishes  for  elastic  spheres.  Even  for  this  simple  mcdel,  the 
complete  density  corrections  are  therefore  not  yet  known. 

For  elastic  spheres!  of  diameter  r0  the  intermolecular  potential  0(r)  is 
given  by: 


4  <  n0 


(6.1) 


It  is  often  convenient  to  consider  0(r)  as  the  limit  of  an  inverse  s-th 
power  law  repulsion: 


(6.1a) 


d2 
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with  r0  x  lim  Xg1/8. 

8-N* 

Bogolubov  has  remarked  that  for  elastic  spheres  one  can  transform  the 
binary  collision  term  Ax(<|,  JjPj.)  of  the  kinetic  equation  to  the  form: 

A,  (If  IF)  =  faj  re  V)  - 

Cf't)>0  (6.2) 

In  here  t  is  an  unit  vector  in  the  direction  of  the  line  of  centers  of  the 
two  spheres  in  the  collision  (P,$i)  ♦($*, ti*)  measured  away  from  molecule  1. 
Equation  (6.2)  has  the  Boltzmann  fo^m,  except  that  the  difference  in  posi¬ 
tion  of  the  two  colliding  molecules  is  taken  into  account.  It  is  the  form 
of  the  collision  integral  from  which  Enskog  starts.*  We  could  therefore 
make  the  transition  to  the  hydrodynamical  stage  similar  to  Chapter  V  but 
starting  now  from  (6.2).  However  we  prefer  to  calculate  directly  the  trans¬ 
port  coefficients  from  the  formula  derived  in  Chapter  V. 

2.  Some  intermediate  results 

Referring  to  the  basic  equation  (5-51)  for  the  determination  of  fx  x, 
we  will  calculate  for  the  elastic  sphere  model  all  terms  in  Lx(p>)  and 
Mlj6>  which  depend  only  on  binary  collisions.  Ore  easily  finis: 

p,  *  -  >  so  p!  *  o 

-  - -f-tf 


•Except  that  Inskog  in  addition  multiplies  the  collision  cross  section  by 
1  4<(5/8)nb#  b  »  (2*/3)  r©*  *  van  der  Waals *  b,  in  order  to  take  roughly  the 
triple  collie lone  into  account. 
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To  calculate  Uj  and  Cjj  we  first  have  to  find  the  integral: 

l4fiejrfA<*‘lr.  ilV 

Calling  this  integral  (Phi*^)*  ve  will  prove  that: 


where  the  notation  is  the  same  as  in  (6.2). 

Proof:  According  to  (4.35b) 

%,i  -  7r« 

Since  the  intermolecular  potential  changes  rapidly  at  r  *  rQ,  there  will  be 
only  a  contribution  to  the  integral  for  r  close  to  rQ.  Therefore  one  can 
neglect  the  dependence  on  px  in  J ri  and  replace  it  simply  by  r^ .  Using 
the  same  cylindrical  coordinates  (I,  b,  0)  for  the  r  integration  as  in 
Section  5»  Chapter  III,  one  then  can  write: 

aZ  .  i  fa,  fa  it  s". 

For  fixed  (b ,0)  S^^6(p-^1  )6(Pi-^a:)  is  a  step  function  of  1  with  discon¬ 
tinuities  at  li  •  ♦  V r0a-b«  and  l2  *  -  -J r02-b2  .  For  I  <  i2 
si^^(P-ni)S(Pi-n2)  *  5(p-Ti  )6(pi-na)  and  I  >  fx  S^Ofp-Tj!  )6(px-qa)  « 
6(p*-Tii  )B(Px*-tja)  where  p*,  px*  are  the  momenta  after  collision.  For 
la  <  I  <  !i  one  can  consider  Sii/5(p-T^1  )B(px-qa)  equal  to  zero  since  the 
momenta  p,  px  are  then  so  high  that  they  will  never  be  close  to  tfc  and  r^. 
In  the  I -integration  the  only  contribution  comes  therefore  from  the  neigh¬ 
borhoods  around  Ix  and  J2.  Let  the  unit  vector  in  the  direction  -r  at 
.  -*• 

(b,  -0,  lx)  be  e,  then  one  sees  that: 
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A;  (1-1,  O  -  -*'*i  ,  Ai  14. 1,1,')-  +*.«; 

i  jftJt  = 

with  (g  •  t)  >  0.  Using  these  results  and  carrying;  out  the  I -integration  one 
obtains  (6.5). 

Introducing  (6.5)  into  the  equations  (5 -kyd,*)  for  %  and  cT^  all 
integrations  can  be  carried  out,  and  one  finds: 


Comparing  this  with  the  general  form  (5.50a,b),  one  concludes  that; 


Next,  from  (5-35 a)  one  can  show  by  a  partial  integration,  that  one  can  write: 


and  therefore  for  elastic  shares  ■  0,  since  -  o. 

Substituting  these  values  in  (5*52)  one  obtains: 

L{(?)  *  -it  +Li(f> 

Mij(f)  •  C#- "  -j£r& if) t 
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where  L*  and  Mjj  are  the  contributions  from  the  triple  collisions  and 
are  given  by: 

Lcf)  =  -  & ** 


(6.7) 


Mg  (?)  =  -  C  rfvjtfA  B„  W. 


Since  the  integral  equations  (5*53  )  are  linear,  one  can  split  the  solutions 
in  a  binary  and  a  ternary  collision  part.  Writing: 


\ll"  =  Vi"*.  +  l/f.w 

U  ’'«•>*“  (6.8) 

,  /«>  .  /')  .  ftn 

Wij  *  Wy  ,  *•.  +  W«^ ,  w 

the  binary  part  will  fulfill  the  equations: 

arc  (»&-*-)&■  (6  9) 

aTX  -  3 

which  except  for  a  constant  have  precisely  the  same  form  as  the  equations 
(3-^3  )  determining  the  zeroth  order  or  Chapman -Enskog  approximation.  There¬ 


fore  one  can  conclude  that: 


Vll 


|/rt) 


,  w;:u  =o 


(6.io) 


As  mentioned  in  the  introduction  the  contribution  of  the  ternary  col¬ 
li  s ions  has  not  been  evaluated. 


Calculation  of  the  t 


rt  coefficients 


For  elastic  spheres  the  Chapman -£ns kog  value  (reference  to  Chapter  10 
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in  Chapman  and  Cowling)  for  the  viscosity  coefficient  is; 

1 1™  =  (i.o/i  ^  (6u) 

According  to  (5*58)  the  first  density  correction  consists  of  two 

parts.  The  second,  or  potential  energy  part  of  can  be  tranpfomed 

for  elastic  spheres  into  the  expression: 


and  from  the  result  (6.3)  for  one  then  finds  that  this  part*  is  equal 

to 

t 

Together  with  the  binary  collision  part  of  V/U),  using  (6.10),  one  obtains 
for  the  total  "binary”  contribution  to  the  density  correction: 


in  agreement  with  Enskog.  The  bulk  viscosity  coefficient  Tjg  vanishes  in 
this  approximation  since  the  trace  of  c^j  is  zero. 

The  Chapman-Enskog  value  for  the  heat  conductivity  coefficient  is  (see 
Chapter  10,  Chapman-Cowling): 

According  to  (5-62),  the  first  density  correction  consists  of  three 
parts.  The  first  potential  energy  part  can  be  transformed  for  elastic 

*It  corresponds  to  the  part  due  to  the  effect  of  "collision  transfer"  in 
the  Enskog  theory.  See  Chapman  and  Cowling,  p.  281  and  p.  282. 
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spheres  into: 


and  hence,  using  (6. 3),  one  finds  that  this  part  is  equal  to: 


Mi i  ^ 


The  second  potential  energy  part  goes  to  zero  for  the  elastic  sphere  model. 
Together  with  the  binary  collision  part  of  losing  (6.10),  one  obtains 
therefore  for  the  total  "binary"  contribution  to  the  density  correction  of 
the  heat  conductivity: 


(6.14) 


again  in  agreement  with  Enskog.  It  should  be  noted  that  the  complete  Enskog 

results  contain  a  rough  estimate  of  the  effect  of  triple  collisions  (see 

* 

footnote  on  p.  83)  which  in  our  terminology  amounts  to  assuming  that: 


W(,)  /  s  < 

wV#)  ~  7 T 


#It  is  again  the  part  due  to  "collision  transfer."  See  Chapman  and  Cowling, 
p.  201  and  p.  287. 


APPENDIX  I 


DERIVAflON  OF  THE  MAYER  EXPRESSION  FOR  j0 


We  have  to  show  that: 


where  the  and  fiz  are  the  qufJititi.es  defined  on  p .  68 . 

pp 

Our  proof  will  be  similar  to  the  one  used  by  Rushbrooke  and  Scoins . 
Using  (5*24)  and  (5.25): 

Rio  -  N +  + . •)]  sv 

Following  Born  and  Green, introduce  a  scaling  variable  /  by: 

l  -If.  ,  V-M' 

then: 

‘fail"* 

Since  the  integrand  vanishes  for  ria  >  r Q/t,  the  dependence  on  I  of  the 
domain  of  the  integration  can  be  neglected.  Therefore 

it  l  (l°  =  firfa  l 

«  PM  C’idl  l  M 

‘  ~3  «•>  *  "  3Pi 
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Next,  one  can  write: 

since  the  integrand  is  independent  of  the  position  of  qx  if  V  is  large.  Us¬ 
ing  again  the  scaling  variable  t,  and  the  symmetry  of  the  integrand: 


-Hh'A:  i  itoiwi*** 

since  again  in  the  q£  and  integration  the  dependence  on  I  of  the  do¬ 
main  of  interaction  can  be  neglected.  Hence  going  back  to  the  unprimed  var¬ 
iables  the  integral  becomes: 

£  &{*-  ttfyi  lll) 

=  -*H*u.u>  -  -tp. 

The  method  can  be  extended  to  the  higher  virial  coefficients,  but  since  they 
are  not  needed  we  omit  the  details  of  the  proof. 


APPENDIX  II 


ISOTROPIC  TENSOR  FIEIDS  AND  NUMERICAL  TENSORS 


Definition  of  an  isotropic  tensor  field 

Let  an  arbitrary  orthogonal  transformation  (proper  or  improper)  of  the 
coordinates  be  denoted  by: 

*  T  4  <*-  Ki  x  Tid^  (!) 

A  tensor  field  is  a  set  of  functions  of  the  which  in  each 

point  transform  as: 

Kij-.-i  ~  TI<T^'*‘  r 

The  tensor  field  is  Isotropic  if  the  transformed  components  are 

the  same  functions  of  the  transformed  coordinates  as  the  original  com¬ 
ponents  were  of  the  original  coordinates,  that  is  if: 


«;x-> 


(?) 


If  (2)  is  only  valid  for  a  proper  orthogonal  transformation,  while  for  a 
reflection  one  has: 


Tie  Tjj  *  T||  K*  •  •  •  y 


one  calls  the  tensor  field  skew-isotropic . 


(2a) 


Definition  of  numerical  tensors 

If  for  an  arbitrary  proper  orthogonal  transformation  T4"  a  tensor  . . ./ 
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satisfies: 

Mij  -l  -  T*  T;l  •••TjJ  tfip—t  (3) 

then  we  call  the  tensor  a  numerical  tensor. 

Note  that  in  the  books  of  Jeffrey2^  our  numerical  tensors  are  called 
isotropic  tensors . 

From  the  definitions  (2)  and  (3)  it  follows  immediately  that  by  inte¬ 
grating  an  isotropic  tensor  field  over  the  whole  space  one  obtains  a  numer¬ 
ical  tensor. 

In  the  following  we  list  some  properties  of  isotropic  tensor  fields 
and  numerical  tensors  which  have  been  used  in  the  text.  For  the  proofs  see 
the  articles  of  Robertson2^  and  the  book  of  Jeffrey.2^®' 

Properties  of  isotropic  tensor  fields 

a)  The  zeroth  rank  or  scalar  field  must  be  a  function  Z(£)  of 

I  *  It!- 

b)  The  first  rank  or  vector  field  must  have  the  form: 

c)  The  second  rank  tensor  field  must  have  the  form: 

Qij(b  -  0,( +  Q,«>S.y 

and,  especially  if  the  trace  of  the  tensor  field  vanishes: 


Properties  of  numerical  tensors 

a)  There  is  no  first  rank  numerical  tensor  except  zero. 
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b)  The  only  second  rank  numerical  tensor  which  does  not  vanish  is  a 
constant  multiplied  by  the  Kronecker  6jj.  Therefore  if  the  trace  of  the 
tensor  vanishes  then  the  tensor  vanishes. 

c)  The  only  numerical  tensor  of  the  third  rank  which  does  not  vanish 

is  a  constant  multiplied  by  e^,  where 

{1  (i,k,|)  is  even  permutation  of  (1,2,3) 

-1  (i,k,|)  is  odd  permutation  of  (1,2,3) 

0  otherwise 

Therefore  if  the  tensor  is  symmetric  in  any  two  indices,  the  tensor  must 
vanish. 

d)  The  general  numerical  tensor  of  the  fourth  rank  is: 

m/l  Siklim  +  /*•  (h  %l«m  +  Sun %kl)  * i  <Sr«i  ” 

with  constants  Especially  if  the  tensor  is  symmetric  in  (i,k)  or  (|,m), 

then 

S«A  Sim  ( Sii  §**  +  Sim  Sm i  -  J  Sit  Sim) 


If  in  addition  the  traces  formed  from  the  first  two  or  last  two  indices  van¬ 
ish,  then 

~  (SuStmi-SlmSn-  ^SikSim) 


Finally  we  will  prove  the  theorem  stated  on  p .  75  •  A  linear  isotropic 
operator  J  acting  on  an  arbitrary  function  f(r|i...^B)  of  the  set  of  vec¬ 
tors  has  the  fonr: 


h>  ■  im-h 


(O 


where  the  kernel  K  has  the  property  that  for  any  orthogonal  transformation 
T: 


kCrplTl-.Tf,)  =  K(?  Il-,p 

Consider  now  the  linear  integral  equation: 

where  Fij...|(p)  is  an  Isotropic  tensor  field  so  that: 

Rij-t  (Tf)  •  TitTjf---Tlr 


Let  us  assume  first  that  the  homogeneous  equation  J(f)  *  0  has  no  solu¬ 
tion  except  f  *  0.  For  each  set  of  values  of  the  indices  i,J,...l  there 
will  then  be  an  unique  solution  of  (6)  which  we  denote  by  f . ,  .  From 

J  •  •  •  I 

(6)  it  follows  that: 

-0Jfr^tK(flV  lr;TV> 

where  we  have  used  the  isotropy  property  (5)  of  the  kernel  K.  From  (6), 

(7)  and  (8)  follows  that: 


-  TiJjf-T„  =  o 


tram  which  one  concludes  that: 
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('Tl-.Tfp 


which  indicates  that  fjj . .  .j  (%  •  •  ^g)  is  8111  isotropic  tensor  field  in  the 
space  of  the  tjj  having  the  same  rank  as  the 

In  the  case  where  ^(f)  *  0  has  solution  we  obtain 

again  an  isotropic  tensor  solution  of  (6)  in  the  function  space  which  is 
perpendicular  to  the  x(*).  Then,  adding  to  this  isotropic  tensor  solution 
the  appropriate  isotropic  tensor  which  is  composed  linearly  of  the  x's,  we 
will  get  the  general  isotropic  tensor  solution  of  (6). 
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